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SUMMARY 


In the first part, the "boundary conditions for an open wind tunnel 
(incompressible flow) are examined with special reference to the effects 
of the closed entrance and exit sections. Basic conditions sire that the 
velocity must "be continuous at the entrance lip and that the velocities 
in the upstream and downstream closed portions must "be equal. For the 
two-dimensional open tunnel, interesting possibilities develop from the 
fact that the pressures on the two free surfaces need not he equal. 

Electrical analogies that might he used for solving the flow in 
open wind tunnels are outlined. Two types are described - one in which 
electrical potential corresponds to velocity potential, and another in 
which electrical potential corresponds to acceleration potential. The 
acceleration- potential analogies are probably experimentally simpler 
than the velocity-potential analogies. 

3n the second part, solutions are derived for four types of two- 
dimensional open tunnels, including one in which the pressures on the 
two free surfaces are not equal. Numerical results are given for every 
case. In general, if the lifting element is more than half the tunnel 
height from the inlet, the boundary effect at the lifting element is 
the same as for an infinitely long open tunnel. 

In the third part is given a general method for calculating the 
boundary effect in an open circular wind tunnel of finite Jet length. 
Numerical results are given for a lifting element concentrated at a 
point on the axis. 


INTRODUCTION 


The basic theory of boundary corrections for an open wind tunnel 
was given by Prandtl many years ago (reference l) and has since been 
used with reasonable success. The infinitely long open jet that was 
assumed in Prandtl 1 s analysis, however, has been frequently questioned 
as an adequate representation for an open wind tunnel, which normally 
has a relatively short jet between closed entrance and exit regions. 
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The present examination of the problem was occasioned by the need for 
boundary corrections for tests in the Langley full-scale tunnel of a 
large helicopter, of which the forward edge of the rotor disk reached 
almost to the mouth of the entrance bell while the rear edge approached 
the exit bell. Previous studies (reference 2) had shown that the 
Prandtl theory was satisfactory for a wing in the usual position in the 
tunnel (about 20 feet downstream of the entrance) j but it was felt that 
this simple theory was inadequate for such far forward and rearward 
locations of the lifting surface, and that some further development was 
desirable . The only previous analysis bearing directly on the problem 
seemed to be that of reference 3, which considered a lifting element 
concentrated at a point on the axis of a circular open tunnel of finite 
jet length J however, the treatment therein was not rigorous, and was 
justified only by a somewhat heuristic discussion, so that its general 
applicability was not obvious. Other studies treated either two- 
dimensional or axially symmetrical conditions (references 4 and 5) and 
also did not consider the closed exit region, so that the extent of 
their applicability to the present problem was not at first apparent. A 
similar German wartime report (reference 6), which did not become available 
until after the present paper was written, would have been more useful 
in this respect because of the generality of its physical discussion. 

Because of the particular shape of the tunnel cross section, a 
reasonably simple solution in terms of available functions seemed unlikelyj 
accordingly, the initial effort was directed toward defining the problem 
in such a way that it could be solved by analogy methods in an electrical 
tank. Identification pf the necessary boundary conditions appeared at 
first to be somewhat perplexing j however, after recognition of some of 
the basic physical phenomena, the boundary conditions were readily 
clarified. The problem is thus considered now to be fairly well under- 
stood, at least insofar as it can be considered linear and uninfluenced 
by turbulent mixing at the free surfaces or by the irregular nature of 
the flow at the exit. As will appear later, however, grave technical 
difficulties exist in the exact solution by electrical-analogy methods, 
so that, for example, actual evaluation of the tunnel interference for 
the large helicopter in the Langley full-scale tunnel, which problem 
instigated the present research, has not yet been accomplished. 

After the boundary conditions were clarified, analytical methods of 
solution were developed for two-dimensional and circular open tunnels. 

These studies have been combined with the discussion of the boundary 
conditions and the electrical analogies to form the present paper, which, 
it is hoped, will serve to clarify basic concepts and establish a sound 
basis for any further work. 

The report is divided into three parts. In part I, the boundary 
conditions are defined and discussed for the open wind tunnel with closed 
entrance and exit sections, and an outline is given of suggested electri- 
cal analogies applicable to the problem. In part II, analytical 
solutions are given for various two-dimensional open- tunnel types. 


NACA TN No. 1826 


3 


together with numerous calculated results. In part III, a method of 
solution for the circular open tunnel is given, together with numerical 
results for the case of a lifting element concentrated at a point on the 
axis of the tunnel. The treatment in every case is a linear one in which 
deformation of the jet houndary is considered to he small. 

The parts were essentially independently prepared. Messrs. Gardner 
and Diesendruck contributed the analysis of part II. Mr. Eisenstadt 
contributed part III. Dr. Katzoff contributed part I, and, in the absence 
of the others, prepared the numerical results of part II, mac^ several 
minor revisions, and served as general editor of the whole. 


I - BOUNDARY CONDITIONS AND 
ELECTRICAL ANALOGIES 


In part I, boundary conditions for an open wind tunnel are discussed 
with special reference to the effects of the closed entrance and exit 
sections. It is shown that the velocity on the free surface is not 
necessarily equal to the velocity far upstream in the closed portion and 
that cross-flows may exist in the free surface, unlike the case of the 
infinitely long open jet. A basic condition - analogous to the Kutta- 
Joukowski condition for the flow at the trailing edge of an airfoil - is 
that the velocity be continuous at the entrance lip. Electrical analogies 
that might be used for solving the flow in open wind tunnels are outlined. 
Two types are described - one in which electrical potential corresponds 
to velocity potential, and another in which electrical potential corre- 
sponds to acceleration potential. 


BOUNDARY CONDITIONS 


Resume of Prandtl's theory .- In Prandtl's original discussion, in 
which the entrance and exit regions are neglected, the tunnel is con- 
sidered as an infinitely long cylinder on the entire surface of which the 
pressure is constant, whence, by Bernoulli's law, the velocity on the 
surface is constant. If this velocity is considered as the sum of the 
undisturbed tunnel velocity U and a small perturbation velocity (u,v,w) 
due to the presence of a body in the Jet, the condition is then 

that (U + u)^ + v^ + w^ » U^ + 2Uu = Constant, from which it is concluded 
that u is constant over the entire surface. Furthermore, since u is 
obviously zero far in front of the body, it must be zero over the entire 
surface . 

A corollary is that, on the jet surface, the perturbation velocity 
normal to u (that is, the circumferential velocity) is also zero, as is 
readily shown from a consideration of the rectangular path SPQR on the 
surface of the jet. (See fig. l(a) .) As has just been shown, the 
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velocity component u parallel to the lines SP and QR is zero} hence , 
the perturbation potentials at points P and Q are the same as at points S 
and R, respectively • If points S and R are far upstream of the body 
their potentials will be equal, so that the potentials at points P and Q 
are equal. The perturbation potential is thus uniform over the entire 
surface, and only perturbation velocities normal to the surface can exist 
at the surface. 

Modification of basic concepts If the closed entrance region is 
near the body* as shown in figure l(b), the preceding discussion and 
conclusions no longer apply. Thus, although u must still be constant 
over the entire free surface, it is no longer necessarily zero} that is, 
the total velocity on the free surface is not necessarily equal to the 
velocity far upstream in the closed portion of the tunnel. The two 
velocities will, in fact, generally be unequal except in special cases 
where equality results from geometrical symmetry of the arrangement. 

(For example, if a horseshoe vortex is located in the horizontal plane of 
symmetry of the tunnel, the values of u at the top and bottom of the 
tunnel would be expected to be equal and opposite: but since u must be 
uniform over the surface, it follows that u = 0.) Furthermore, the 
velocities in the jet surface normal to u (that is, the circumferential 
velocities) are, in general, no longer zero (except for axially symmetri- 
cal flows, such as that produced by a source on the axis of a circular 
tunnel) so that two surface points at the same longitudinal position, 
as P' and Q 1 (fig. l(b)) do not necessarily have the same values of the 
perturbation potential. 

Entrance-lip condition .- Consider, for simplicity, the symmetrical 
case of figure l(b) , in which the lifting element is on the horizontal 
plane of symmetry of the tunnel. Since u is zero on the free boundary, 
the perturbation potential is constant along the elements AB, CD, 

EF, • • ., although, as just indicated, it is not necessarily the same 
for all these elements. This one boundary condition for tjie open section - 
that the potential be constant along each of these elements - does not 
suffice, however, to define the problem uniquely. In fact, as will be 
obvious from the subsequent discussion of electrical analogies, the 
potentials of these elements may be quite arbitrarily assigned without 
violating this condition or the boundary condition on the closed portion 
of the tunnel (that the normal derivative of the potential be zero at the 
wall) . In order to avoid this lack of uniqueness, further conditions 
must be sought. The most important of these is that the velocity be con- 
tinuous (in particular, not infinite) at the entrance lip (points A, C, 

E, . . .) . This condition takes cognizance of the fact that, because of 
viscosity, the physical flow leaves the lip smoothly, just as it leaves 
the trailing edge of an airfoil} the condition is, in fact, strictly 
analogous to the Kutta- Joukowski condition for the trailing-edge of an 
airfoil, which similarly takes into account the basic viscosity effect 
and provides uniqueness where otherwise an infinity of solutions would exist. 
It is recognized that, just as the Kutta- Joukowski condition does not 
always suffice to predict airfoil lift very accurately, the corresponding 
condition for the open tunnel may similarly oversimplify the entrance- lip 
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flow; however, as with the airfoil, the condition is probably adequate 
where the flow is not subject to an excessive pressure rise on approaching 
the lip. Keferences 4, and 6 used the condition, and reference 6, in 
addition, discussed it from the physical viewpoint and compared it with 
the airfoil trailing-edge condition. 

Concerning the downstream end of the open section, the exit lip may 
be considered to correspond to the leading edge of an airfoil and no 
effort need be made in an idealized flow analysis to eliminate infinite 
values of u at this edge. 

Jet contraction or expansion .- It has already been pointed out that, 
with a body in the jet, the velocity on the free surface is not necessarily 
equal to the velocity far upstream in the closed portion. During the 
course of the investigation, it was noted that solutions could be obtained 
showing a difference between these two velocities, even when there was no 
body in the jet. Such a flow corresponds merely to a contraction or 
expansion of the jet, as indicated in figure 2. Thus, in figure 2(a), the 
velocity on the free surface is lower than the upstream velocity and 
remains so even as it approaches the exit, in spite of the gradual contrac- 
tion of the jet, because of the continuously increasing surface curvature. 
The velocity suddenly increases at the exit lip and finally is established 
at a value greater than that of the upstream velocity. With reasonable 
ratios of entrance to exit area, the flows of figure 2 may be readily 
obtained experimentally. 

The significance of this expanding or contracting flow is that it 
represents a solution that satisfies all the boundary conditions previously 
discussed and is nevertheless undesirable. In order to avoid such 
solutions, a further condition must accordingly be recognized; namely, 
that the velocities in the closed portions far upstream and far downstream 
of the open section be equal. 

It may be objected that in the normal design of an open wind tunnel 
the exit section is made larger than the entrance section. The purpose 
of the increased area is to allow for the reduced velocity toward the 
surface of the jet resulting from turbulent mixing with the surrounding 
still air. Increasing the exit area by other than the correct amount will 
result in the type of flow indicated in figures 2(a) or 2(b), with a 
corresponding velocity gradient along the center of the tunnel. In any 
potential-flow solution these viscous effects cannot be considered. 

Spillage .- When an airfoil is tested at a high lift coefficient in 
an open tunnel, the downward deflection of the jet may result in appreci- 
able spillage from the lower lip of the exit, together with lack of contact 
of the main flow with the upper lip. (See fig. 3(a).) The air lost by 
spillage is replaced by air (of, however, a lower total pressure) 
entrained in the exit. Even without otherwise considering the distortion 
of the free surface, these flow characteristics might seem too much at 
variance with the previously assumed characteristics to permit applica- 
tion of the theories being discussed. The calculations of part II for 
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the two-dimensional open tunnel (that is, a rectangular tunnel with 
closed sides but open top and bottom) show, however, very little differ- 
ence between the tunnel- indue ed-downwash distributions for the tunnel 
with two exit lips and the tunnel with one exit lip. That is, if 
figure 3(a) is assumed to represent a two-dimensional flow, the fact 
that the upper lip of the exit is out of the flow field so that the 
lower lip takes over the entire burden of straightening the jet does not 
greatly affect the induced downwash. 

The effect of the exit lip on the flow phenomena is the least clear 
of the various phases of the present problem. For open wind tunnels 
having essentially unflared exits, similar to that Indicated in figure 3(a), 
the suggestions of the preceding paragraph are probably adequate. The 
exit of the Langley full-scale tunnel, however, has a large bell mouth, 
and when airplanes are being tested at high lift coefficients a downward 
deflection of the air off the lower part of the bell, roughly as indicated 
in figure 3(t) , occurs. Whether the previously suggested concepts or, 
indeed, any linear theory can serve satisfactorily for this case seems 
questionable . 

Unequal surface pressures .- An interesting method of avoiding 
spillage suggests itself in the case of the two-dimensional open tunnel: 

If the space below the tunnel is inclosed, an excess pressure will be 
built up in this space, compared with the pressure in the space above 
the upper free surface, so that the flow will be pushed up sufficiently 
to eliminate the spillage and ensure precise contact of the lower free 
surface with the lower exit lip. (See fig. 3(c).) The extent to which 
a free two-dimensional jet can be deformed by a pressure difference 
across its boundaries, or, stated differently, the extent to which a two- 
dimensional free jet will deform in order to follow the only available 
path, is indicated by the smoke-flow photograph in figure 4. The setup 
consisted merely of a two-dimensional open jet with entrance and exit 
sections displaced vertically relative to each other, arranged between 
transparent side walls, and provided with enclosed spaces above and below. 

Details of interest in the figure, in addition to the jet 
deformation, are: 

(1) Separation of the flow from the upper lip of the exit, because 
of the large angle of entry. A small bell mouth at the exit lip might 
have prevented such separation. 

(2) The rough flow on the upper surface at the entrance, compared 
with the smooth flow on the lower surface, reflects the fact that the 
boundary layer approaching the entrance is subjected to a rising pressure 
on the upper surface and a dropping pressure on the lower surface . 

(3) Because of turbulent mixing at the free surfaces, a certain 
amount of the air in the closed chambers above and below the jet is 
entrained in the jet. An equivalent quantity must be released, or 
skimmed off, at the exit in order that the total quantity in each chamber 
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remain constant. This circulating mechanism results In the apparent 
overflow at the two exit lips. The return of the skimmed- off part to 
the jet surface can he seen at the bottom of the photograph. 

• 

Tunnel without a closed exit .- Some mention is made in the subse- 
quent discussion of the hypothetical open tunnel having a closed upstream 
entrance region hut no closed exit region, the open section thus extending 
downstream to infinity. Calculations for such an arrangement (see part II) 
are generally simpler than for the actual tunnel with the closed exit, 
and give very nearly the same answer, provided that the region of 
interest is much closer to the entrance than to the exit, as is usually 
the case. For this arrangement, solutions with an arbitrary contraction 
or expansion of the jet cannot exist, so that no effort need he made to avoid 
them. The solution for the general unsymmetrical case, however, will 
show the jet velocity downstream at infinity to he different from the 
velocity upstream, in the closed part. The possibility that, in the two- 
dimensional case, different pressures might he assumed on the two free 
surfaces still exists for this type of tunnel, hut the resulting jet will 
have a constant curvature after leaving the neighborhood of the body. 

An upstream condition for the "infinitely long” open tunnel and a 
correction to the results of reference 8 .- In many discussions of the 
two-dimensional open tunnel, the set of images indicated in figure 5(a) 
is used to satisfy the boundary condition that u = 0, and the resulting 
flow shows an upflow in front and an equal downflow in hack, with no 
induced downwash at the wing itself. Actually, however, if the jet 
issues from a horizontal closed entrance - no matter how far upstream - 
it will remain essentially horizontal (because it is not subjected to 
any vertical force) until it reaches the wing. (See fig. 5(t>)») In 
order to eliminate the undesired upstream upwash, a uniform downwash 
should therefore be added to the solution indicated in figure 5(a) • 

(compare reference 7 , p. 30^*) Addition of this downwash does not 
affect the boundary conditions, since u is still zero at the boundary. 

This case is discussed quantitatively in part II, where it is shown that 
the entrance-lip condition automatically provides the correct answer. 

Among the rectangular wind tunnels for which corrections were 
given in reference 8 is a type with closed sides but open top and bottom. 

The calculated corrections for approximately square cross sections are 
approximately equal to those for the completely closed tunnel, a sur- 
prising result in view of the absence of any top or bottom constraint. 

The result is actually in error, as was discovered in an experimental 
effort to verify it (reference 9) » In seeking to explain the errors 
the author of reference 9 pointed out that the image system used in 
reference 8 should have included an infinite row of vortices at infinity, 
and he showed how, by taking into account this row of vortices, the 
correct, answer could be obtained. It could not be shown, however, that 
the extent of this row of vortices is of a higher order of infinity than 
is their distance from the origin, as is necessary if their effect is to 
be considered. The method of the preceding paragraph thus appears to be 
much simpler and more rigorous in such cases than is a discussion of the 
image vortices at infinity. One simply observes that the image system 
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of reference 8 provides an angle correction factor 5 of 0.25 for the 
flow far upstream of the wing , whereas 8 should he zero far upstream; 
a correction of -0.25 should therefore he added to all values of 6 
computed by this image system for points within the tunnel. 

Su m m ar y of boundary conditions .- A basic physical characteristic 
of the flow is provided hy the condition that the velocity he continuous 
at the entrance lip, which also helps to provide uniqueness. The 
velocity on the free surface is not necessarily the velocity far upstream 
in the closed portion; in fact, for the two-dimensional case, it is even 
possible for the pressures on the two free surfaces to he different from 
each other. Equality of the velocities in the upstream and downstream 
closed portions has been recognized as an additional condition. 

Neglecting the upper portion of the closed exit may he desirable if the 
flow is so depressed that it does not make contact with the upper part 
of the exit. Neglecting the entire closed exit region may appreciably 
simplify the problem without introducing excessive inaccuracy if the 
region of interest is much closer to the entrance than to the exit. 

In general, adequate treatment of the exit (for large lift on the body 
in the tunnel) seems very unlikely. 

The discussion in the preceding sections has concerned mainly the 
physical flow conditions, and relatively little interpretation in terms 
of boundary conditions on the perturbation potential has been given, 
although such formal interpretation would appear a trivial task. The 
reason that this extension has not been made is that, in a number of 
instances, as will appear subsequently, slight modifications of the basic 
viewpoint, leading to somewhat modified boundary conditions, are desirable 
for convenience of solution. Accordingly, the statements of the boundary 
conditions on the perturbation potentials will be given when the solutions 
are discussed. 


SUGGESTED ELECTRICAL ANALOGIES 
Velocity- Potential Analogies 


Basic concepts of the analogies .- In the analogies to be discussed 
in. the present section (none of which have yet been constructed) , the 
perturbation velocity potential in the space within the wind tunnel is 
considered analogous to the electrical potential in a dilute electrolyte 
solution contained in a vessel of the same shape. An insulating material 
such as Bakelite , the conductivity of which is negligible compared with 
that of the solution, provides a boundary where the normal potential 

gradient ^ is zero; and a metal, the conductivity of which is practi- 
cally infinite relative to that of the solution, serves as a constant- 

potential boundary along which the longitudinal gradient ^ is zero. 

ox 

In such a setup, current is analogous to velocity except for a difference 
in sign (in the usual convention, current flows down a voltage gradient 
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whereas air flows up a velocity-potential gradient ) } in order to remove 
this difficulty, the sign convention for electrical potential is reversed 
in the following discussion. 

For greater clarity of exposition, the two-dimensional analogies are 
treated in detail, the three-dimensional analogies appearing as reasonably 
obvious extentions or modifications. It will be remembered, however, that 
any application will be found in three-dimensional problems, inasmuch as 
most of the two-dimensional problems can be solved analytically. 

A two-dimensional vortex may be represented by two long metal plates 
separated by a thin insulator. (See fig. 6(a).) A flow corresponding 
to a vortex located at the edge of the plates is set up by applying a 
difference of potential across the plates. If the perturbation flow 
that results from the presence of the vortex in the tunnel has a hori- 
zontal velocity component, this representation is no longer adequate 
because it requires the potential to be uniform along each plate. Rigor 
in this case would require that the plates be composed of a number of 
separate sections, with each pair separately activated. (See fig. 6(b).) 

In this way it is possible to provide a potential difference between 
upper and lower surfaces that is everywhere equal to the desired circula- 
tion, without requiring that the potential be uniform along the entire 
upper surface or lower surface. A horizontal velocity component normally 
occurs only when the lifting vortex is asymmetrically located in the 
tunnel. For simplification, only the simpler representation of figure 6(a) 
is used in the remaining sketches . 

The element of lift in three-dimensional flow is the horseshoe 
vortex of zero span, which is the same as a semi- infinite line of doublets. 
It may be represented by a pair of long narrow metal strips separated by 
an insulator. (See fig. 6(c) •) As in the two-dimensional analogy, if 
the lifting element is asymmetrically located in the field the strips 
must be made up of short pieces, with each pair separately activated. 

The horseshoe vortex of finite span is represented as in figure 6(d), 
provided there are no appreciable perturbation velocities in its plane. 

Evaluation of interference velocities .- The vertical velocity 
component in the tunnel corresponds to the vertical voltage gradient in 
the electrolyte, which can be determined by measuring the voltage differ- 
ence between a pair of short wire electrodes mounted one above the other 
a fixed distance apart. The tunnel interference at any point is found by 
measuring this voltage difference (relative to that across the two plates 
representing the vortex) first in the simulated tunnel and then in a 
large tank for which the boundary interference is either negligible or 
so small that it can be adequately computed by simple methods. Since 
the theoretical flow field for the second case is known, the ratio of 
these two gradients, together with the distance from the pair of wires 
to the lifting vortex, should suffice to evaluate the boundary interfer- 
ence. Tlie distance between the pair of wires need not be measured 
because only the ratio of the gradients is required. Similarly, the 
exact design and dimensions of the simulated lifting vortex are of no 
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significance, provided the gradients are determined at reasonable 
distances from it. In general, "boundary interference at the vortex itself 
cannot be found directly by this method, but may be determined by 
interpolation between or extrapolation from neighboring points. 


Two- dimensional closed-open tunnel .- For simplification of the 
nomenclature, the open tunnel with closed upstream region but without a 
closed exit is designated the closed- open tunnel. The open tunnel with 
closed upstream and downstream regions is designated the closed- open- 
closed tunnel. 


Figure j(a) illustrates the setup for a two-dimensional closed-open 
tunnel with a vortex on its center line. Shaded lines indicate 
insulating boundaries, where ^ = 0, and heavy unshaded lines indicate 

metal boundaries on which $ is constant. The upstream closed portion 
should be so long that the potential is essentially uniform at its 
upstream end} the length indicated on the figure should suffice. The 
open region should similarly be so long that the vertical flow between 
the vortex strips and the boundary strips no longer changes with distance 
downstream} again, the length indicated on the figure should suffice. 
From the condition of velocity continuity at the entrance lips and the 

fact that is zero along the free boundaries, it follows that ^ 

ox ox 

must be zero at the edges of the two closed boundaries. The potentials 
on the two free boundaries must, therefore, be adjusted until the differ- 
ence between the potential of each and the potential of a thin feeler 
electrode Just upstream of its edge is zero. For the symmetrical condi- 
tion shown, the single variable voltage source indicated will provide 

at both edges simultaneously . 


zero 


Figure 7(b) illustrates the setup for the two-dimensional closed- 
open tunnel with the vortex in an off-center position. A single variable 
voltage source across the two free boundaries is now no longer capable 
of simultaneously satisfying the continuity condition at both edges, so 
that an additional variable voltage source and an upstream electrode 
are required. The current in the closed part of the tunnel flowing into 
this upstream electrode corresponds to an upstream perturbation velocity. 
This upstream perturbation velocity constitutes the previously mentioned 
difference between the velocity far upstream in the closed part and the 
velocity on the free surface. The concept here is slightly at variance 
with previous discussion, which considered a perturbation velocity along 
the free surface, with the far upstream velocity appearing as the 
undisturbed velocity U. As the analogy is set up, however, no pertur- 
bation velocity may appear* along the free surfaces because they are at 
constant potential} hence, the total velocity on the free surfaces must 
be considered as the undisturbed velocity U and any difference between 
this velocity and the velocity far upstream appears as an upstream 
perturbation velocity. As appears in part H, this viewpoint is also 
found convenient in the analytical solution of these problems. 
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Two- dimensional closed-open- closed tunnel .- The setup for the two- 
dimensional closed-open-closed tunnel with a vortex on the center line 
(fig. 8(a)) is an obvious modification of the corresponding setup for 
the closed-open tunnel. The same would he true for the off-center vortex 
except for the necessity of satisfying the condition that the velocities 
in the upstream and downstream closed regions he equal. Thus, the setup 
of figure 8(h) provides an upstream perturbation velocity hut no down- 
stream perturbation velocity, and cannot, therefore, solve the problem 
completely. An additional flow, found by the setup of figure 8(c) must 
be included. An electrode is here located at both the upstream and 
downstream ends, and the potentials relative to the free boundary are so 
adjusted that the entrance-lip condition is satisfied. It is apparent 
that in order to satisfy this condition the downstream current flow will 
be much greater than the upstream current flowj that is, the downstream 
perturbation velocity for a contracting or expanding jet is much greater 
than the upstream perturbation velocity. Because of this difference, a 
suitable amount of the flow of figure 8(c) may be added to that of 
figure 8(b) to produce equal upstream and downstream perturbation 
velocities. 


Displacement of the free surfaces .- The current density normal to 
the surface of a metal plate representing a free surface is proportional 

to ^ — and corresponds to the local vertical perturbation velocity. 

The total vertical displacement at a point on the free surface is then 


given by ^ dx integrated from the entrance lip to the point. In 

J W 

particular, the integral along the entire lower free surface of a closed- 
open-closed tunnel represents the displacement at the exit lip and it 
may be measured by means of an ammeter in the line that goes to the 
lower metal plate . 


If the pressure on the lower free surface can adjust itself so 
that the displacement at the downstream end is zero, the perturbation 
velocity at the lower surface will be different from that at the upper 
surface. If the perturbation velocity on the upper surface is taken as 
zero, that on the lower surface will be negative, so that the potential 
on the lower surface must drop uniformly from entrance to exit. Such a 
variation could be accomplished if the lower surface were represented by 
a number of short metal strips instead of a single plate. 

Three-dimensional closed-open and closed- open- closed tunnels . - The 
analogies for the three-dimensional tunnels are obvious modifications of 
those for the two-dimensional tunnels . The free boundary may not be 
simulated by a single metal cylinder because, as was previously noted, 
different elements of the free boundary do not have the same potential, 
although the potential is constant along each element. The free boundary 
must thus be simulated by a number of longitudinal metal strips, 
insulated from each other, with a feeler electrode immediately in front 
of each. When the lifting element lies in the horizontal plane of 
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symmetry, the entrance-lip condition may he satisfied without an additional 
electrode . If the lifting element is not in the horizontal plane of 
symmetry, an upstream, electrode will he needed, with the potential of 
each strip adjusted relative to this electrode. For the unsymmetrical 
closed- open- closed analogy, the requirement that upstream and downstream 
velocities he equal necessitates further measurements with a setup 
corresponding to that of figure 8(c). 


Acceleration- Potential Analogies 

Basic concepts of the analogies .- The pressure has the properties 
of a potential - designated acceleration potential - in a field consisting 
of a small perturbation flow superposed on a uniform stream. If the 
pressure in the undisturbed stream Is taken as zero, then the perturbation 
velocities are related to the pressure by the following equations: 


u 







i_ r x dp 

puj.oo dy 


JL_ P x dp 

P u J „ S 

— 00 


dx 


dx 


_P_ 

PU 


where 

p density 

t time 

p pressure 

Since, by the first equation, u is proportional to p, it is simpler 
merely to consider the perturbation velocity u itself as the potential, 
with v and w given by the following equations: 


and 



8u 


dx 


w 



du 


dx 
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The necessity of performing an integration in order to determine v or w 
is a "basic disadvantage of the acceleration-potential analogy compared 
with the velocity-potential analogy in which v and w are measurable 
directly « 


In the analogies to he discussed in the present section, the pertur- 
bation velocity u is considered analogous to the electrical potential 
in a dilute electrolyte solution. A metal serves as a boundary along 
which u is constant, and the local intensity of current flowing into 

it gives an insulator serves as a boundary where is zero. The 

on dn 

lifting element in either two or three dimensions is represented by a 

pair of short metal plates separated by an insulator; when the upper 

plate is maintained at a higher potential than the lower plate, the 

arrangement represents a thin airfoil with suction (large u) on its 

upper surface and pressure (small or negative u) on its lower surface. 

The current at each of the two plates should be the same in order that 

the slope of the airfoil surface (proportional to v) be the same on 

both upper and lower sides. In order always to satisfy this condition 

the voltage source activating the lifting element should not be tapped 

to any other electrode in the field. 


Two-dimensional closed-open tunnel .- The setup for the two- 
dimensional closed- open tunnel with the lifting element on the center 
line is shown in figure 9(a) . The walls of the upstream closed region 

are represented by insulators, which establish that ^2=0 at every 

r x n °y 

point; hence, the condition that v = I 

J-00 oy 


dx = 0 at every point on 


the closed boundary is satisfied. The two free boundaries are repre- 
sented by metal; and electrically connecting them, as shown, satisfies 
the further condition that they have the same potential (the same u) . 
The flow of current into the lower boundary then equals the flow of 
current out of the upper boundary, so that the ultimate downstream 
value of v will be the same on both upper and lower boundaries, as is 
desired. In fact, for the symmetrical case illustrated, the value 
of v will be the same at all pairs of opposite points on the two free 
boundaries; so the boundaries will be everywhere parallel. No special 
attention need be paid to the entrance lips - the entrance lip condi- 
tion is automatically satisfied since the potential u is continuous 
at these points (although the potential gradients at these points are 
infinite) . 


For the off-center position of the lifting element (fig. 9(b)) no 
modification of the circuits is needed. The difference between the 
potential in the upstream closed region and the potential of the free 
boundaries, which is the upstream perturbation velocity u, is measured 
with the aid of the probe P. As in the symmetrical case, the ultimate 
downstream value of v will be the same for both the upper and the 
lower boundaries; however, it is no longer true that the two boundaries 
will be everywhere parallel, and the ultimate width of the jet will be 
different from the width of the closed part. 


Ik 
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Two- dimensional closed- open- closed tunnel .- The value of v at the 
downstream, end of the free boundary is given by the total flow of current 
into the metal plate and, in general, is not zero. At the lip of the 
closed exit, however, v must be suddenly reduced to zero in order for 
the flow to follow the solid boundary ) hence, a short electrode must be 
added at the exit lip, and as much current must be forced out of it as 
flows into the long electrode that represents the open boundary} that is, 

the integral of ^ dx along the free boundary must be canceled at the 
oy 

exit lip. The setup (fig. 10(a)) therefore shows a voltage source to 
supply this current and means for measuring and equalizing the current 
flow into adjacent electrodes . If these additional short electrodes are 
omitted, the setup will correspond to a tunnel the exit section of which 
has been alined with the deflected Jet (fig. 10(b)) because the .condition 

that ^ = 0 on the closed exit boundary would merely permit v to 
oy 

remain at the value it had at the end of the free boundary. 


For the off-center position of the lifting surface, a similar setup 
is used and, as before, probes in the regions far upstream and far down- 
stream are used to determine the potential u in these regions relative 
to the potential of the free boundary. Since these potentials far up- 
stream and far downstream will not be equal, an additional perturbation 
field must be provided such that the sum of the two fields will have the 
same potential in the two regions. This additional perturbation field, 
which corresponds to a contraction or expansion of the Jet, is provided 
by the setup shown in figure 10 (c) . It is clear from this figure that 
the downstream perturbation potential is much greater than the upstream 
perturbation potential} this result corresponds to that indicated in the 
velocity- potential analogy. 


The condition in which the pressure on the lower free surface is 
higher than that on the upper free surface is easily represented by 
applying a voltage difference between the two surfaces. (See fig. 10(d).) 
The corresponding displacement of the lower surface, however, is not so 


readily obtained. The vertical velocity at every point is y ~ dx, so 

n a>u 


that the displacement at each point is 


¥ 


dx dx. 


rau 
J dy 

In order to accom- 


plish this integration ^ must be determined at points along the 

oy 

boundary, perhaps by breaking the long plate into a number of short 
pieces and determining the current flowing into each. 


Three-dimensional closed-open and closed-open-closed tunnels .- The 
analogies for the three-dimensional tunnels are again obvious modifica- 
tions of those -for the two-dimensional tunnels. For the closed-open 
analogy, the free boundary may be represented by a single cylinder of 
metal (fig. 11(a)). For the closed- open- closed analogy, the free boundary 
must be represented by a number of separate strips (fig. 11(b)) in order 
that the total current into each strip may be measured and an equal 
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current forced out of the short strip Immediately "behind it. Contraction 
or expansion of the jet is represented as in the two-dimensional casej 
hut the setup that would correspond to different pressures along differ- 
ent strips seems to have no practical significance in the three- 
dimensional case. 

Some form of this acceleration- potential analogy is probably the 
most convenient for solving problems similar to that of the helicopter in 
the Langley full-scale tunnel. Simply neglecting the exit, as with a 
closed- open tunnel, permits the free surface to be represented by a 
single sheet of metal and eliminates any measurements of current flow to 
or from the surface. Improved accuracy should be attainable by cutting 
the sheet into two parts with two short strips at the rear. (See 
fig. 11(c).) The need for many strips seems unlikely, at least in view 
of the previously mentioned uncertain definition of the physical flow 
in the region of the exit. 

Correspondence between velocity- potential and acceleration-potential 
analogies As has already been indicated, the acceleration potential is 
identical with the x-component of the perturbation velocity and is hence 
merely the x-derivative of the perturbation- velocity potential. It is 
of interest to point out the related fact that the acceleration potential 
analogies are, in a sense, the x- derivatives of the velocity- potential 
analogies. For example (see fig. 12), 

(1) For the velocity-potential analogy, an infinitely long double 
layer represented a lifting element located at its forward edge. The 
difference between two such double layers, of which one is shifted 
slightly relative to the other, is merely the short double layer that 
was used in the acceleration- potential analogy. 

(2) For the velocity- potential analogy, the free boundary consisted 
of constant-potential strips on which the potentials were so adjusted 
that the gradient was zero at the leading edge. If each strip is now 
shifted and subtracted, there remains a long strip, with a short strip 
at the front and back. Since, in the velocity- potential analogy, the 
gradient was zero at the entrance lip, the short strip at the front may 
be neglected. The remainder corresponds to the arrangement used in the 
acceleration- potential, analogy, and the fact that the total current after 
the subtraction must be zero corresponds to the fact that the total 
current out of the short strip must be made equal to Hie total current 
into the long strip. 

( 3 ) When the lifting element was off-center, the velocity- potential 
analogies required electrodes upstream and downstream, with uniform 
current flow along the upstream and downstream closed regions . That the 
subtraction eliminates these current flows corresponds to the fact that 
no upstream or downstream electrodes are used in the acceleration- 
potential analogies. 
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Technical Difficulties 

It should he pointed out that the analogies here described may he 
rather unwieldy, experimentally. Sven for the simplest types cf analogies 
the literature indicates considerable uncertainty as to the most satis- 
factory electrolyte and electrode materials, and appreciable difficulty 
in balancing capacitances (alternating current is generally used in 
analogies, in order to minimize polarization at the electrodes) . In the 
present analogies, the need for separate current sources that are exactly 
in phase and the large capacitances that will certainly characterize the 
vortex and the open-boundary representations should greatly complicate 
the technique. Perhaps the use of direct current instead of alternating 
current, with nonpolarizing electrodes (as platinized platinum), would 
be a more practical approach in this respect. Simultaneously satisfying 
the entrance-lip condition at a number of points around the inlet (or 
satisfying the corresponding exit condition for the acceleration- 
potential analogies) may also turn out to be very difficult. 


RESUME OP PAET I 


The most significant points of the preceding discussion of open 
wind tunnels and their electrical analogies are as follows: 

1. Continuity of velocity at the lip of the entrance cone is a 
basic characteristic of the flow in an open wind tunnel. 

2. Equality of the velocities in the upstream and downstream closed 
regions is a further condition on the tunnel flow if extraneous longi- 
tudinal pressure gradients are to be avoided. 

3* The velocity on the free surface need not be the same as the 
velocity in the closed upstream region. In general, the two velocities 
are the same only when the lifting element lies in the plane of 
symmetry of the tunnel. 

4. For the two-dimensional open tunnel the velocities on the two 
free surfaces need not be equal. If the space below the lower free 
surface is closed off, the pressure on the lower free surface will adjust 
itself so that the displacement at the exit lip is zero. 

5 . Considerable uncertainty exists with regard to conditions at the 
exit or the mathematical equivalents of these conditions. Correspondingly 
certain compromises in complying with the idealized downstream conditions 
may be justified in a determination of boundary interference. 

6 . In any analysis that neglects the closed entrance and exit regions 
the condition of zero upstream induced flow must be retained. 
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7 • Electrical analogies of either the velocity-potential or the 
acceleration- potential ‘type may he devised to correspond to most of the 
problems discussed. 

8. In electrical analogies that represent velocity potential, by 
electrical potential, the condition of continuity at the entrance lip 
appears troublesome, especially for three-dimensional tunnels j however, 
the exit conditions are easily represented. 

9* In electrical analogies that represent acceleration potential 
by electrical potential, the entrance- lip condition is automatically 
satisfied, but fulfillment of exit conditions is troublesome. Rough 
approximation of the exit conditions may, however, be adequate for many 
purposes . 

10. Acceleration-potential analogies are experimentally simpler 
than velocity-potential analogies. 


II - TWO-DIMENSIONAL TUNNELS 


In part II, boundary- induced velocities in two-dimensional open 
tunnels are derived with special reference to the effects of the closed 
entrance and exit regions . The cases treated are : 

(1) Tunnel with a closed entrance (upstream) region but without a 
closed exit region 

(2) Tunnel with a closed entrance region but with only one exit 
lip (corresponding to a condition in which the downward deflection of 
the flow is so large that the flow makes contact only with the lower 
exit lip) 

(3) Tunnel with closed entrance and exit regions 

(4) Same as case 3, "but with different pressures on the two free 
surfaces 

Numerical results are given for all cases. 


SYMBOLS AND DIMENSIONS 


Each tunnel is idealized as a strip of uniform height h, having a 
stream velocity V, and containing a point vortex of strength r. For 
simplification of the present development, lengths and velocities will 
be made nondimens ional by dividing by h and V, respectively, and the 
vortex strength will be made nondimens ional by dividing by hV. 
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Essentially, then, the solutions will be developed for a vortex of 

strength r' = in a tunnel of unit height} and in the following 

list of symbols the lengths in the complex planes are in terms of h, 
and the complex velocities are in terms of V: 


h 

tunnel height 

V 

tunnel velocity 

i 

complex variable of physical plane (5 + ii)) 

h 

location of vortex in | -plane 

z 

complex variable of transformed plane (x + iy) 

Z 1 

location of vortex in z-plane 

q 

complex velocity in physical plane (u - iv) 

Q 

complex velocity in transformed plane (u - iv) 

A, B, C, M, N 

real constants 

r 

vortex strength 

r* 

nondimensional vortex strength 

qi(^!) 

induced complex velocity at £ when vortex is 

siUi) 

induced complex velocity at £-j_ 

a 

abscissa of exit lip in transformed space 

W 1 

a complex velocity in the foim of an elliptic : 
of the first kind 

w 2 

a complex velocity in the form of an elliptic 
the second kind 

W 3 

a complex velocity 

l 

variable of integration 

G 

function defined by equation (2) 
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K, K 

complete elliptic integrals of the first kind, with 
modulus l/a 

E, E' 

when not followed by parenthesis, complete elliptic 
integrals of the second kind, with modulus l/a j 
with upper limit indicated in parentheses, 
incomplete elliptic integrals of the 
second kind, with modulus l/a 

F, F' 

incomplete elliptic integrals of the first kind, with 
modulus l/a, and with upper limit indicated in 
parentheses 

R.P. 

real part 

I.P. 

Imaginary part 

c 

airfoil chord 

c l 

airfoil lift coefficient 

€ 

tunnel- induced angle, radians 

*b 

horizontal perturbation velocity at free boundary 


Sub script 


i 

induced 


BOUNDARY CONDITIONS 


The two-dimensional tunnels discussed are considered to have their 
fixed and free "boundaries parallel to the real axis, with the main tunnel 
flow from, left to right. The physical plane (in which lengths and 
velocities have "been made nondimensional as just described) will be 
designated the £-plane, with the complex perturbation velocity u - iv, 
or q(£), subject to the following conditions: 

(1) on each fixed (or closed) boundary, I.P.q(£) = -v = 0 

(2) On each free (or open) boundary, R.P.q(^) = u = 0 or a constant 

(3) At each lip of the closed entrance section, q(£) is continuous 
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CASE 1 - CLOSED-OPEN TUNNEL 


Total perturbation velocity * - By the transformation 

z = e n ^ (l) 

the tunnel in the plane, represented By an infinitely long strip of 
unit height, is transformed to the upper half of the z-plane . The 
correspondence between points is shown in figure 13 • 

The complex velocity (rather than the more usual complex potential) 
is considered to be retained in the transformation, and the problem is 
thus to find a function Q(z), where 

u - iv = Q(z) = q( t) 


such that 

(1) On the closed sections of the boundary, that is, for z real 
and | z | < 1, 

I.P.Q(z) = 0 

(2) On the open sections of the boundary, that is, for z real 
and |z| > 1, 

E.P.Q(z) = 0 

(3) For z = +1, 

Q(z) = 0 

(4) Q(z) is finite at infinity 

Consider the complex velocity G(z) corresponding to a vortex at z^_ 
and its reflection at z"is 



This function, which is of order l/z^ at infinity, satisfies co ndi- 
tions (l) and (4) but not (2) and (3) • Functions of the form \/l - z 2 , 

z \1 - z^ , . . . satisfy conditions (l), (2), and (3)« At infinity, 

\Jl - z2 is of order z and z \|l - z^ is of order z^j therefore. 
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either of the products G l/l - z^, Gz\/l - z^, or a linear combination 
of the two satisfies the four conditions and has a pole of the first 
order at z-j_. These facts suggest that the desired velocity function 
is of the form. 


Q(z) = i(-U - 1 _\ a + Bz) \fl - z 2 (3) 

\z - z-l z - z-y) 

where A and B are as yet undetermined real constants. 

The values of A and B are to be determined such that the pole 
at Z]_ represents a vortex of strength - T ' in the plane . Thus 

r 1 =j£q(0 d£ =$Q(z) dz =i^Q(z) 

where the integral is taken about the point zj_. By Cauchy's integral 
formula 

r' = -2 (A + Bzj. / 1 ~ Zl2 00 

Z 1 

The values of f and z-y are known, so that this complex equation 
can be solved for the two real constants A and B. Substituting 
these values in equation (3) will thus give the desired complex velocity 
function. 

Tunnel- interference velocity.- The tunnel- interference velocity is 
defined as the difference between the total perturbation velocity q(£) 
due to the presence of the vortex in the tunnel and the velocity due to 
a vortex in an unbounded medium. That is, the tunnel- interference 
velocity giU^l) is 


<3i(Ui) = (A + Bz) N/ 1 ' z2 + 2*u i - $!) (5) 

u t 1 

If the vortex is on the axis of the tunnel, that is, if + 7^, 

then from equation (l), = iy^, and equation (4) gives 
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A = 0 


2\/l + yi 2 

If the point of evaluation is also on the axis (that is, £ = 4 + 
whence z = iy), then 


( 6 ) 



\ 


y 2 + l 
71 2 + 1 


+ 


ir' 



(7) 


Thus, if the vortex is on the axis, the interference velocity at all 
points on the axis has only a vertical component. 

The interference velocity at the vortex itself is the limit of 
expression (5) as z approaches z]_. The term containing z - zp 
offers no difficulties and its limit is readily evaluated ’ 


lim ; (a + Bz)\/l - z 2 = - 7 j~“(a + Bz-i)\/l - zi 2 = t : 

z^z 1 z - z x d ?i 4y]_ 


where the last equality follows from equation (4) . The remainder of 
equation (5), after Jt(£ - £1) is replaced "by log — , is 


(A + Bz)\nT- z 2 + 
z - zi 


2 log — 
zi 


(A + Bz)\jl - z 2 log -r- + z - z-l) 


= i 


zi 


(z - Zl )log ^ 


This expression is of the form for z = zj_. Differentiating numerator 

and denominator, according to L 'Hospital's rule still leaves both equal 
to zero at z = zi (that the derivative of the denominator is zero 
at z = Z]_ is obvious ) that the derivative of the numerator is also 
zero at z = zp can be verified with the aid of equation (4)). A second 
differentiation yields the following expression for the limit 

as z approaches z^_: 


NACA TN No, 1826 


23 


-(a ♦ ♦ 2BZ1 ^ 


/rr 


Z 1 


zi c 


2z n 


This fraction can he greatly simplified hy use of equation ( 4 ), and the 
result, added to the previously derived limit, gives the desired correc- 
tion at the vortex: 


lim q ± (g,5l) = qjftl) = i 

5^5i 


r_|_ r 1 

k + 2(1 - Zi 2) 



+ 


r'zi 

*nyi 


(8) 


For the special case in which the vortex is on the tunnel axis (z = iy) , 
this expression reduces to the following form (after substituting for B 
from equation (6)): 


< 3 i 



yr 


1 + yr 


(9) 


Upstream perturbation velocity If the vortex is not on the tunnel 
axis, A will not be zero (compare equation ( 6 )). The tunnel interfer- 
ence velocity far upstream in the closed part of the tunnel is found by 
putting £ = - 00 and z = 0 in equation (5), which then reduces to 




2Ayi 


which is real. For this unsymmetrical case, therefore, a finite longi- 
tudinal perturbation velocity is found far upstream in the closed part. 
As was pointed out in part I, such results appear because the problem 
was set up so that the longitudinal perturbation velocity on the open 
boundary is zero. If the velocity far upstream in the closed part is 
to be taken as the base, the result means merely that the longitudinal 

velocity on the open boundary exceeds this base velocity by ^yi . 

I z i! 2 

that there is a corresponding difference in pressure between the closed 
part and the space surrounding the jet. 
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Limiting case of completely open tunnel.- With increasing distance 
of the vortex from the closed entrance (that is, with increasing yp) , 
expression (9) approaches - • As was pointed out in part I, the 

image system normally used to satisfy the boundary condition on an 
infinitely long, open, two-dimensional tunnel produces no induced flow 
at the vortex itself, and only after introduction of the additional , 
condition that the upstream flow he horizontal is this value of - 

for the induced-flow correction obtained. In the present development, 
however, it is seen that the condition of continuity at the entrance 
lips automatically takes care of this condition on the upstream flow 
direction, even when the entrance and the vortex are infinitely far 
apart. 

No further discussion of the completely open tunnel will be given 
here inasmuch as this case has been adequately treated by the method of 
images. (See reference 6, p. 302.) 


CASE 2 - TONNEL WITH ONE FIXED EXIT BOUNDARY 


Perturbation velocity.- As before, the transformation z = e rt ^ 
transforms the tunnel, considered as an infinite strip of unit height, 
into the upper half of the z-plane. The correspondence between points 
is shown in figure l4. The conditions on the complex velocity Q(z) are: 

(1) On the real axis, I.P.Q(z) = 0 for -1 <z <1 and for z >a 

(2) On the real axis, E.P.Q(z) = 0 for z <-l and for 1 < z < & 

(3) For z = ±1, Q(z) = 0 

(4) Q(z) is finite everywhere in the upper half-plane except 

at z = a and at z = zp. As noted in part I, Q(z) will be infinite 
at z = a 

'The function 0 (z), given by equation (2), will again be used as a 
factor that is real along the entire real axis and has the desired type 

of singularity at zp. The functions 811(1 Z V' ^ ~ ~ chan 8 e 

from pure real to pure Imaginary as z passes through ±1 or through a, 

and, furthermore, are zero at z = ±1 and infinite at z = a. At 

infinity they are of order z 1 / 2 and z^/ 2 , respectively. Theref ore, 
as before, a linear combination of G(z)^^— 3 and G(z)zy^ _ ^ 
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satisfies the preceding conditions and has a pole of the first order 
at zp. Q(z) is therefore of the form: 


Q(z) 




( 10 ) 


The real constants A and B are determined by the same condition as 
"before, which gives 


r' 


2 (A + Bz-^ 
Z 1 




( 11 ) 


Tunnel-interference velocity .- The interference velocity qi(C^i) 
is the difference Between the perturbation velocity and the velocity due 
to a vortex located at the same point in an unbounded medium: 


<3i(£,£l) 



7"bi)< A + 


q - z 2 ir 1 
a - z 2n(t - £-]_) 


( 12 ) 


This expression does not simplify appreciably if the vortex is located 
on the axis, and the interference velocity at a point on the axis due 
to a vortex on the axis is not normal to the axis . 

The interference velocity at the vortex itself is the limit of 
expression (12) as z approaches zp. Proceeding as in the preceding 
case gives: 

r 1 r'a 1 1 - zi 2 r'zj 

2(l - Zi 2 ) 4(a - zi) V a - zx 4iyi 

It may be shown with the aid of equation (ll) that this equation reduces 
to that for the closed- open tunnel as a goes to infinity. 


qi(£i) = 1 


CASE 3 - CLOSED- OPEN-CLOSED TUNNEL 


Perturbation velocity .- The transformation z = e n £ transforms the 
tunnel space into the upper half of the z-plane with correspondence 
between points as indicated in figure 15. The boundary conditions on 
the complex velocity Q(z) are: 
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(1) On the real axis, I.P.Q(z) = 0 for |z| <1 and for | z| >a 

(2) On the real axis, R.P.Q(z) = 0 for 1 < |z| < a 

(3) Por z = ±1, Q(z) = 0 

(4) Q(z) is finite everywhere in the upper half-plane except 
at z = +a and at z = zp 

(5) Q(0) = Q( 00), (this equation corresponding to the condition 
noted in part I that the perturbation velocities in the upstream and 
downstream closed regions he the same) 


The function G(z) given by equation (2) is again used as a factor 


of Q(z) . The functions 


- *2 


I a 2 - z 2 ' 


'a 2 - z 2 ' 


and z c 


- z 2 


' 2 2 
a - z 


satisfy 


Q 

conditions (l) , (2), and (3), and are of orders 1 , z, and z at infinity, 
respectively. By the same reasoning as before. 


Q(z) = i 


£ - zp z - zp 


^t-Va + Bz + Cz 2 ) 


a 2 - z 2 


(Ik) 


The condition that the pole of the first order at zp represents a 
vortex of strength I” is: 


r’ = - 


2 (A + Bzp + Czp 2 ) 1 - zi 2 


2 2 
a^ - zp 


(15) 


Condition (5) is satisfied by equating the two foims of equation (l4) 
for z equal to zero and equal to infinity. Thus 


Q(0) 



. ~ g yi a 

l z i|^ a 


lim Q(z) 

Z - x» 


lim — 
z^~ ( z 


-2yi 

zp) ( z - zp) 


(A + Bz + Cz 2 ) 



= -2ypC 
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whence, by condition (5) 


- g yi a 

z i | 2 a 


-2yiC 


or 


C 



(16) 


The complex equation (15) and the real equation (l6) suffice for 
evaluating the three real constants A, B, and C in equation (l4) • 

Tunnel- interference velocity.- The tunnel interference velocity is 


- 1(7^ 



(A + Bz + Cz 2 ) 


1 - z 2 ir 1 

- z^ 2 it ( £ - e 


(17) 


If the vortex is on the tunnel axis, that is, if Z]_ = iyi, equation (15) 
gives 


B 



C 





J 


(18) 


Comparing this last equation with equation (l6) , which reduces to the 
following form for Zj = iy^ 


C 



shows that for this symmetrical case 


A = C = 0 


28 


NACA TN No. 182 6 


since, as is clear from figure 15, a ^ 1. The interference velocity at 
a point on the axis (z = iy) due to a vortex on the axis is thus 


<3i 




ir'yjy l + y 2 a 2 + y x 2 ir* 

y 2 - y^ 2 V a 2 + y 2 1 + y^ 2 2jt ( 5 Si) 


(19) 


which is normal to the axis. 

The interference velocity at the vortex itself is the limit of 
expression (17) as z approaches z Proceeding as before gives 



r' + r'*i s 

4 2 


a' 


2 _ 


(a 2 - Zl 2 ^l - z x 2 ) 


+ (B + 2Cz!) 


V 2 - 


Z 1 


r»zi 

(20) 


For the case in which the vortex is on the axis (zp = iyjJ the 
normal velocity at the free boundary ( z = x, where 1 < x < a) is given 

iy 


or 


Q(x) = ii 


iyi x + iy 


Bx 


L - x 2 
a 2 - x 2 


Q(x) 


2iyi 
x 2 + y x 2 


•Bx 



JL_ 

X 2 


CASE k - CLOSED-OPEN- CLOSED TUNNEL WITH UNEQUAL 
PRESSURES ON THE FREE SURFACES 


( 21 ) 


Boundary conditions.- As indicated in part I, the two-dimensional 
closed- open- closed tunnel may develop unequal pressures on the two free 
surfaces if a closed space exists below the lower free surface. Within 
the limits of the present linear theory, this pressure difference corre- 
sponds to superposing on the flow discussed in the preceding section an 
additional perturbation velocity field Q(z) that 

(1) Has no singularities within the tunnel 

(2) Satisfies the condition of continuity at the inlet lips 
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(3) Has a horizontal component equal to, say, +1 on the upper free 
boundary and -1 on the lower free boundary 

(4) Has no vertical component on the closed boundaries 

(5) Is zero at infinity upstream and downstream 

Rewritten as conditions on the complex velocity Q(z) in the z-plane 
(fig« 15 ), these conditions become 

(1) Q(z) has no singularities in the upper half of the z-plane 

(2) Q(z) =0 at z = ±1 

(3) On the real axis, Q(z) = -1 for 1 < z < a, and Q(z) = +1 
for -a < z < -1 

(4) On the real axis, I.P.Q(z) = 0 for |z | < 1 and for | z | > a 

(5) Q(z) = 0 for z = 0 and for z = 00 

Outline of method*- Consider the following two functions of zi 



They can be considered as complex velocities having the following 
properties along the real axis (compare reference 10) : 

W 1 is real between 1 and -lj between 1 and a, or between -1 and -a, 
its real part is constant but an imaginary part is introduced^ beyond a 
or “ a , the imaginary part is constant while the real part approaches 
zero} R.P.w^z) «= -R.P.w^-z); I-P-v^z) = I.P. Wl (-z) 

w 2 has the same properties as w^_ except that beyond a an ^ -a 
its real part approaches 00 and - 00 , respectively 
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Maps of the two functions are shown in figure 16. It should 
obviously he possible to find a linear combination of these two func- 
tions, Mw]_ + Nwg, such that for 1 < z < a 


R.P.(mw]_ + Nwg^ = ~1 


for -a < z < -1 


and beyond a or -a 


R.P.(Mwp + Nwg) = +1 


I .P . (Mw]_ + Nwg^ = 0 

A simple additional function wo to be discussed subsequently is needed 
to satisfy the condition at infinity. The desired velocity function for 
the closed- open-closed tunnel with unequal pressures is thus of the form: 

Q(z) = Mw^(z) + Nw 2 (z) + w^(z) 

The constants M and N are derived in the two following sections . 

Evaluation of integrals .- In the following development, the modulus 
of all the elliptic integrals is l/a j the modulus will therefore not be 
indicated in the symbols. In the designations for the incomplete 
elliptic integrals, E, E 1 , F, and F 1 , the terms in parentheses are the 
upper limits of integration. Then 


wi(l) = R.P.w^(a) = 


dz 


\|(! - z 2 )(a 2 - z 2 ) 


= ^K 
2 \ a 


wgCl) = R.P.w^Ca) = 




I.P.w-^(a) = 


dz 


ii v(i - - * 2 ) 


which by the substitution z^ = a^ - (a^ - 1 jl^ reduces to 



dZ 


1 - ! 2 ¥l - 


iK' 
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I.P.WgCa) 




dz 


p P 

which by the substitution z = a 




reduces to 




= iaK 1 - iaE 1 


Solution of simultaneous equations for M and N .- With the aid of 
the four formulas Just derived, the two previously mentioned equations 
in M and N may be written 


% + NaE = -1 
a 


+ NaK' - NaE' =0 
a 

which are easily solved simultaneously for M and N. By introducing 
the following relation between the complete elliptic integrals (refer- 
ence 11, p. 520) 


EK* - KK' + KE' = | 

the expressions for M and N are finally obtained in the following 
forms : 


M = — (K 1 - E') 
Jt v 


N = - 


2K 1 

ait 
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Value of Mwi + Nw^ at infinity.- The constants M and N have 

"been determined so that I.P.(Mwi + Nw 2 ) = 0 at infinity $ furthermore, 
R.P.Mwi =0 at infinity, as is clear from figure l 6 . Therefore, the 
value of Mwi + Nu 2 at infinity is merely R.P.Nv^ at infinity. It 
is necessary to investigate this limit before choosing the form of wo, 
because, as was previously noted, the purpose of wo is to 
provide Q(z) =0 at infinity. The limit may be written 


R.P.Nw^C 00 ) 



dz + 



dz 


The first term is simply NaE. In order to evaluate the second term, 
substitute z = 7 

L 




The first term is -NaK. In order to evaluate the second term, it is 
noted (reference 12 ) that 


d_ \ fi 1 - Z 2 )(a 2 - Z 2 ) 
dZ 7 


Z 2 n/(i - Z 2 )(a 2 - f) 

+ + 

Z 2 V(l - Z 2 )(a 2 - Z 2 j 7(1 - Z 2 )(a 2 - Z 2 ) \J(l - Z 2 ) (a 2 - Z 2 ) 
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Transposing terms in this equation gives 


a2 

l 2 \J (l - l 2 ) (a 2 - Z 2 ) 


a2 

\f C 1 ■ z2 X &2 ■ z2 j 





whence 


Na 2 f 1 - dZ - 

Jo - 1*1*2 - l*) 



dZ 

V (i - * 2 X &2 - ?2 ^ 



IT \/( 1 - Z 2 Xa 2 - Z 2 ) 
Z 


1 

0 


= NaK - NaE - N- 



1 

0 


The first two terms on the right are exactly canceled ty the two terms 
previously obtained, so that the final result is: 


lim E.P.NwgCz) 
z — >oo 


lim 

1-pO 



1 

0 


= lim 

Z — >00 



= lim Nz 

Z —* x > 
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Derivation of w^.- In order to cancel the effect of the 

terms Mw]_ + Nw£ at infinity, the function W3 must approach -R.P.Nwg (z) 
as z increases without limit. In addition, it must have no singularities 
in the upper half of the z-plane, it must satisfy the condition of con- 
tinuity at z = ±1, it must he a pure real on the fixed boundaries and a 
pure imaginary on the free boundaries , and it must be zero at z = 0 . 

It is readily formulated as 


w 3 


(z) 



That this function satisfies the first condition is readily shown by 
writing it in a slightly different foim and expanding the radical: 


lim -Nz 



1 


lim -Nz 
z— > 00 




= lim -Nz(l + — ~ 1 + . . . 

z— > 00 ^ 2z 


N(a 2 - 1 ) 

= lim -Nz 1 — - L -... 

2z 

z — *• 00 

Comparison of this expression with that for lim R.P.Nw^Cz) shows that 

z — > 00 

the difference between the two expressions approaches zero as z 
approaches infinity. That the function satisfies the remaining conditions 
is readily verified by inspection. 

The complex velocity function for the closed- open- closed tunnel with 
unequal pressures is, finally. 


Q(z) 



(K' 



z 


dz 

- 7% 1 - T) 



( 22 ) 
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or 


Q(z) = |(K' - E')F(z) - ^pE(z) 


2K' |l - z 2 

a 2 * Z tTI. 

a 2 


where the modulus of the elliptic integrals is l/a. 


(23) 


Induced velocity on the axis .- For the special case in which z = iy, 
the preceding equations for Q reduce to a somewhat simpler form. The 
procedure will "be only outlined here, inasmuch as the manipulative steps 
are similar to those already described. 


Replacing 
tuting y 2 = 


z with iy in the expression for 

o 

- a reduces the first term to 


wi 


and then substi- 



E 1 ) 


K' 



The same substitutions, together with the previously described technique 
from reference 12 reduces the second term to 




2iK , < 

x_ / 1 + y 2 

7t 

W 1 + £ 


_ * a 2 


K' 




The third term is found directly as 


m. * JUjE 

a2 « V 1 + 

1 a2 

The total simplifies to the form 



( 24 ) 
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Normal velocity and normal displacement at the free surface .- The 
normal velocity on the free surface may he written in the following form: 


x 


I.P.Q(x) 

l<x<a 


= ii(E* - E') 

it 


dx 


Jl 


f ■ 1 X 1 - i 



x 2 

1 " „2 2iK' 


x 2 


x 2 - 1 


x 2 

1 - ~ 


By the substitution of x 2 = a 2 - (a 2 - l) 2 2 , the integrals are readily 
reduced to standard forms of incomplete elliptic integrals, a n d the 
equation takes the form 


I.P.Q(x) 

l<x<a 


21 





k'e 



(25) 


The normal displacement, or distortion, of the free surface is 
found by integrating this expression along the free surface in the 
physical plane: 


Normal displacement at x = j I.P.Q(x) d| 


= / I.P.Q(x) ^ 

Jl 

The integral may be evaluated numerically j however, the third term 
of I.P.Q(x) is amenable to analytical treatment: 


_ 2iK' 



2iK' 




E' - E 


i a 


- x 2 


Va 2 - 1 / 


K' - F ,: 


a 2 - x 2 

a 2 - 1 
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At the edge of the exit lip, where x = a, this expression reduces to 



NUMERICAL RESULTS 


In the following sections are described some numerical results that 
were computed by the preceding equations in order to show the magnitudes 
of the entrance and exit effects. It will be noted that, since the 
complex velocity has been made nondimens ional by dividing by V, the 
component v is identical with the tunnel- induced angle e, in radians, 
and the component u is the fractional increase in the horizontal 
velocity. The equivalence of the two ordinate scales indicated in the 
plots of the results follows from the equation 

I = 

Y 2 

Closed-open tunnel. - In figure 17 are shown calculated values of 
the induced downwash angle along the tunnel axis for various positions 
of the lifting vortex along the axis. The figure shows that for = 1.0 
and 1*5, the induced angles at the f vortex itself (I = 1.0 and 1*5, 

respectively) are almost exactly which is the value for an infinitely 

long open tunnel j and, furthermore, the two curves are symmetrical about 
the point I = ip. In fact, within the accuracy of the plot, these two 
curves are identical with the curve for an infinitely long open tunnel. 

It may be concluded that the closed entrance has no effect if the vortex 
is more than one tunnel height from the entrance. For = 0 .5 which 
is a more likely location of the wing, the induced velocity at £ = ^ 
is 0.48r', and the curve is no longer exactly symmetrical about the 
point | = £]_; however, these differences from the conditions for the 
infinitely long open tunnel are too small to be practically significant, 
so that the usual inf ini te- open- tunnel theory is still adequate 
for = 0.5* For less than 0 . 5 , the deviations from infinite- 

open-tunnel theory become larger rapidly, until, when the vortex is in 
the plane <jf the entrance = 0Y the induced angle at the point £ = £, 
is only 

k similar discussion applies for the vortex in the closed portion 
of the tunnel < oj, although this case normally has no practical 
significance. For £]_ = -1, the induced angles in the neighborhood of 
the vortex are practically identical with those for an infinitely long 
closed tunnel; however, in the open region (I >0), the curve is consider- 
ably different from that for the infinitely long closed tunnel (shown as 
the dashed curve in fig. 17 ) • 
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Symmetrical closed- open- closed tunnel .- In figure 18 are shown 
similar curves for a symmetrical closed- open- closed tunnel of which the 
length of the open section is 1.5 times the tunnel height. All curves 
show a sharp reduction in the induced angle as the closed exit is approached 
and entered^ however, for Ip < 1, the closed exit has practically no 
effect on the induced velocities at and upstream of the vortex. The 
induced angle at the vortex decreases rapidly as the vortex moves down- 

r 1 

stream from about 5p = 1.0, and is -jj- in the plane of the 
exit (£p = 1.5) . 

Closed- open- closed tunnel with one exit lip .- Figure 19 shows 
results for a tunnel similar to that just discussed except that one exit lip 
is omitted. The two curves shown are very similar to the corresponding 
curves for the symmetrical condition. As was pointed out earlier, the 
horizontal component of the induced velocity on the axis is not zero for 
this unsymmetrical configuration. Values of this horizontal component 
have been plotted in figure 20 for the same two vortex locations as in 
figure 19* The values are seen to be very small in the forward part of 
the tunnel but become quite large in the neighborhood of the exit lip. 

The effect is consistent with the concept of the exit lip as a concen- 
tration of vortices having total strength equal and opposite to that of 
the bound vortex and serving thereby to turn the air back to its original 
direction. The fact that the two curves are practically identical lends 
further support to this viewpoint. 

Comparison of the three tunnel types .- In figure 21 are compared 
the induced- angle curves for Ip = 0.5 and 1.0 for the three tunnel 
types just discussed. It is seen that the differences are slight up to 
about’ £ = l.Oj beyond this value the curves for the closed-open tunnel 
continue to rise, while the others descend rapidly. The effect of the 
closed exit is somewhat larger for the tunnel with two exit lips than 
for the tunnel with one exit lip. Although the induced angles become 
slightly negative in the downstream closed region they eventually return 
to zero. 

Symmetrical closed- open- closed tunnel with unequal pressures on the 
two free surfaces. - By means of equation ( 2 k) calculations were made of 
the induced vertical velocities on the axis of a closed-open-closed 
tunnel of jet length equal to 1.5 times the tunnel height and having 
equal and opposite horizontal perturbation velocities (-ub and ub) on 
the upper and lower free surfaces, respectively. The results are plotted 
in figure 22. The curve shows that the vertical velocity component (or 
the induced angle) has an almost linear variation along the axis, which 
corresponds to the fairly uniform curvature of the jet that would be 
expected to result from the pressure difference between the upper and 
lower surfaces • For this same condition, the integral of the normal 
velocity along the free surface (equation (25)), which is the downward 
displacement of the jet boundary at the exit lip, was found to 
be 3 «89ub . 
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For a vortex r' located at £l = 0.5 in the equal- pressure case, 
the integral of the normal velocity along the free surface (equation (21)), 
which is proportional to the upward displacement of the jet boundary at the 
exit lip, was found to be 1.20r'<> Accordingly, zero displacement at the 
exit, corresponding to the existence of a closed space above or below the 
jet (reference 2), will result if the flow described in the preceding 
paragraph is superposed on the equal- pressure flow in such proportion 

that 3.89ub = 1.20r'j that is, = 0.31* The corresponding effect on 
the induced angle at £ = li is found as follows: at I = li = 0.5, 


for the equal- pressure case (fig. 18) is 0.48. From figure 22, — at 

£ = O.5 for the unequal-pressure case is -1.44. Since 0.31 X -1.44 = -0.45, 
it is seen that, if spillage at the exit lip is prevented, the induced 
velocity in the region of the vortex is nearly eliminated. A similar 
comparison of the slopes of the curves in figures 18 and 22 in the 
neighborhood of I = 0.5 shows that the induced curvature in the region 
of the vortex is also nearly eliminated. 


Resume of numerical results .- The induced angle at the lifting 
vortex is essentially that for an infinite open jet if the vortex is more 
than half the tunnel height from the entrance and the exit. The induced 
angles for case 2 (one fixed exit boundary) are nearly the same as for 
case 3 (symmetrical exit), so that any failure of the flow to contact the 
upper exit lip should not appreciably affect the tunnel correction. 
Finally, for case 3 , if enough of the different- pressure flow is added 
to assure zero displacement of the free boundary at the exit (that is, 
if spillage at the exit is prevented, as by enclosing the space into 
which the spillage would normally occur) , the induced angle at the 
vortex may be nearly eliminated. 


Ill - CIRCULAR TUNNELS 


In part III an outline is given of a general method for calculating 
the boundary effect in an open circular tunnel of finite Jet length. 

The solution, involving expansions in Bessel functions, is somewhat 
similar to the solution for a closed circular tunnel (reference 13) , 
but is constructed so that it satisfies the condition of uniformity of 
pressure over the open boundary and also the condition of continuity of 
velocity at the entrance lip. Numerical results are given for a lifting 
element on the tunnel axis. 
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SYMBOLS 


ft, i, 5 


ft, p, e 


a, b 


3 

q 

$ 

3>n 


m 


u 


( 2 ) / , 
8m (ft) 


^ran 


(J) 


^sm 

(l) 


L m 


( 2 ) 


m 


rectangular coordinates in units of the tunnel radius 
with origin at lifting element (see fig. 23 ) 

cylindrical coordinates (see fig. 23 ) 

I -coordinate of entrance and exit lips, respectively 

variable of Integration 

variable of integration (see reference 13) 

disturbance potential associated with body (or with 
vortex system) 

tunnel- Induced potential 

tunnel- induced potential in closed circular tunnel 
residual potential (<J> - 

Bessel function of the first kind of order m 

constant longitudinal perturbation velocity on free surface 


g m ^(|) m^ Fourier sine coefficient of 


Sp" 


P=1 


"fcll 

m Fourier cosine coefficient of 


4“V> 

n coefficient in series for 


8m 


a *A 

3p” 

(J) U) 


P =i 


s th zero of J m ' (not including the zero at the origin) 


m.^ 1 Fourier sine coefficient of 




+ § 


dt 


’£ 1 


p=i 


m th Fourier cosine coefficient of - ^(^9 




P=1 


tunnel- induced velocity parameter 




NACA TN No. 1826 


kl 


p density of fluid 

L lift of lifting element 

w tunnel- induced velocity normal to the 111 -plane 

V free -stream velocity 

R tunnel radius 


ANALYSIS 

Introduction 


In the analysis of the three-dimensional, circular, closed-open- 
closed tunnel, an appreciable simplification results when the tunnel axis 
lies in the plane of the horseshoe vortex . For off-center locations of 
the horseshoe vortex, or for a source- sink body on the axis, or for the 
general unsymmetrical disturbance, certain complications arise that are 
related to the fact that the pressure on the free boundary is then not 
equal to the pressure at ±°° in the closed parts of the tunnel. That 
is, for these cases, if the net perturbation velocity is zero far 
upstream and downstream in the closed parts of the tunnel, a constant 
longitudinal perturbation velocity u ^ 0 will exist on the free surface. 
(See parts I and II.) A similar complication results for a source in a 
completely closed tunnel. 

The analysis described in the following section is applicable 
directly to the case in which the tunnel axis lies in the plane of the 
horseshoe vortex and for which the longitudinal perturbation velocity on 
the free surface is zero. (See part I.) In the succeeding section are 
derived the additional terms needed for the solution of the more general 
problem. The significance of the titles of these two sections will 
become clear in the analysis. 


Cylindrically Symmetric Term Omitted 

Boundary conditions and formal egression for <J>a ,_ The solution 

i s _ developed in cylindrical coordinates (|,p,9) where the i- axis 
coincides with the tunnel axis and 9 is measured from the horizontal 
plane . The relations of these coordinates to the rectangular coordi- 
nates (i,T|,t) are indicated in figure 23- The distance variables |, q, £ 
and p are considered in units of the tunnel radius. 

Let $q(|,p,0) be the disturbance velocity potential associated 
with the lifting body in unlimited space (in particular, the velocity 
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potential of a horseshoe vortex) . It is desired to find a func- 
tion O(t } p,0), harmonic inside the cylinder p = 1, for which (see 
part I) 


8(3>q + $) 


dp 


= 0 


p=l 


(4 < a, £ > h) 


d($o + $) 
81 


P=l 


= o 


(a < | < b) 


where the region i < a, I > b is the closed portion of the tunnel 

and a < I <b is the open portion of the tunnel. In addition, according 

to the condition of continuity at the entrance lip (part I ) , the 

must be continuous at I = a. The function $ is 

P=1 

then the velocity potential of the additional flow due to the tunnel 
boundary . 


derivative Si 
dp 


The function $ is conveniently considered in two parts: 


$ = + 0 A 


where Oq is the known tunnel- induced potential for the same vortex 
system in a completely closed circular tunnel (reference 13) • The deter- 
mining conditions for $ A are then 


( 1 ) 

( 2 ) 

(3) 

(4) 


A3> a = 0 


d*A 

ST 


= 0 


P=1 


yr 


p=i 


d(frp + $C) 
8| 


P=1 


8o a 

8T 


= 0 


P=1 


(p < 1) 

(6 < a, I > b) 
(a < £ < b) 


(6 = a) 
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Hie function $a H^y 1© expressed formally (see reference 13) as 


$A = 


sin m0 — 


1 r r* ( 1), * 

«J iiJ m '(iq) dq / 8m l ; (P) cos q(P - |) dp 


m = 1 


n l f 

+ cos m 9 - — , , , 


(26) 


where 


fc (1, <o 


and 


(0 


are, respectively, the m 


th 


cosine coefficients of the Fourier series for 


d$ A 

dp 


and 


sine and 
q and P 


P-1 

are variables of integration. The integrations over p would, in 
general, have the range to +«>$ however, condition (2) shows that 

the functions g^J)(p), j = 1, 2, are zero from -00 to a and from b 
to +°°. The convergence of this function and its derivatives to the 
desired function $a and its derivatives is discussed in the appendix 
of reference 13. A modification is necessary because of the disconti- 

0<ft A 

nuity in that may exist at | = b . For this case , the desired 

ap p=l 

convergence may be proved for regions bounded away from the 
circle p = 1, i = b. 


The assumption of zero perturbation velocity on the surface of the 

jet is equivalent to the assumption that the expansion of - — (^Q *’7^') 

o * 


in a Fourier series in 0 contains no term independent of 0. For this 
reason no m = 0 term appears in expansion (26) . The next section 
discusses the somewhat special treatment that is required when the 
Fourier series contains a term independent of 0. 


P=1 


Evaluation of gn7^ ( £) .- The function $a given in the preceding 
equation satisfies conditions (l) and (2) regardless of the precise form 
of the functions g m ^^(^). It is now desired to find the func- 
tions for a < 5 < b such that % will satisfy conditions (3) 

and (4) . To this end the functions are represented by infinite series 
of the form 


= h mO (j) 3in f b^i 



hj^(j) sin niti-=-S 
b - a 


(27) 
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Since gm^ (a) thereby equals zero, condition (4) is automatically- 

satisfied o Values of j) (t>) are here assumed to he finite, instead 
of infinite (see part I)j the corresponding inaccuracy, however, is 
considered to exist mainly in the immediate region of the exit 
lip (|=b,p=l). 

Substituting this series in equation (26) gives 


_.£» 


$ A = 


b ffln ^"*"^Pmn(|,p) sin m0 + hmn ^ ^ P^p ( E , o ) cos m0 


m = 1 n = 0 


( 28 ) 


where 


v>a,p) -r 


1 Jm(ip<3) 


— / 1 — t 1 7 • \ d-q ? f - - a cos q(p - |) dp 

n J 0 iqJ m '(iq) J a 2 b - a s/ p 


and, for n ^ 0, 


r ,, , 1 / " JmdP?) , r b , B 

Wi,p)-;7 o Tijyrrq) d<! J a 3ln 


- a 


h - a 


cos q(p - I) dp 


In the evaluation of these two expressions, the inner integrals may be 
found analytically and the outer integrals, which converge rapidly, may 
be found numerically. It is possible, however, by means of contour 
integration similar to that discussed in reference 14, to transform the 
infinite integrals into infinite series that are more convenient for the 
present purpose. The contour integration and the resulting infinite 
series are given in appendix A. 

Differentiating equation (28) with respect to I and taking p = 1 

gives 


SF 


P=1 


%00 00 



sin m0 


^mn^A) 


+ ^mn 


( 2 ) 


cos m0 


^mn( £>1) 

8! 


(29) 
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The constants h^-n ^ are to he determined so that condition ( 3 ) is 


satisfied" For this purpose the function 
in a Fourier series in 6 : 


_ d ($ 0 + $c) 


ST 


is expanded 


'P=l 


- 


0 ^ + $, 


S I 


3£l 


p=i 


m = 1 


r m ^(|) sin m0 + r m ^ 2 ^(^) cos m 9 


(30) 


Equating coefficients in equations (29) and (30) in order to satisfy 
condition ( 3 ) then gives 


% (j) (d = 2^ v (J) %r^ (j - 1, 2) 


n = 0 


It is assumed that the functions *m^(l) can he satisfactorily approxi- 
mated hy a finite number of terms of these series. This assumption seems 

reasonable, inasmuch as — 5 B j ^ ■ ’ — is hounded as n approaches infinity 


(see appendix B) and hum 0^ approaches zero as n approaches infinity. 
Thus . 


r n (j) (i) * XI V 0) 

n = 0 


The functions r m ^(0 and T j. are compu -5ed at a set of 

points (ijj, i = 0, 1, 2, . . . I, where Iil, The coefficients 
are then determined (method of least squares) so that the expression 


1=0 


N 


r m <J>(?i) - XI V (J) 

n = 0 0 


is a minimum for all values of m and J. For each pair of values of m 
and j, this condition gives N + 1 equations for the N + 1 
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unknowns 


h m0^> h ml^ 


timN 


(j) 


as the N + 1 partial derivatives 


with respect to l%n^ must he equal to zero. These equations 


are 



i = 0 


r n (J) Ui) v (J) 


n = 0 


dT 


-W 

(t = 0, 1, 2, 


= 0 


. . N) 

Remarks on the computations The points ft -A and the value of N 
are chosen so that the addition of more points "and increasing the value 
of N will no longer appreciably affect the results. It is clear that 
the point I = h cannot he used and that care must he taken not to 
choose too large a proportion of the points ft^ in the neighborhood 
of t = a and I = bj any such attempt to describe more accurately 

the infinite values of g^^U) at | = h or of its derivatives 

at I = a and I = h with a finite number of coefficients will cause a 

large error in the approximations to g m (I) elsewhere in the inter- 
val a < | < b. 

Since the functions r m ^ (| ) rapidly approach zero as m approaches 
infinity, the preceding equations need be solved for only a small number 
of values of m. The values of hmn ^ thus obtained can be used to 
give an approximation to the function $a (equation (28)): 


M 


N 


*A = 


m = 1 n = 0 


1 %n^ p mn(£,p) sin m0 + itaa^ 2 ^%n( t,p) cos mflj 


(31) 


Any desired interference velocity may now be obtained by differentiating 
this series term by term and adding the results to the corresponding 
interference velocity for the closed tunnel. 

The vertical induced velocity in the plane of symmetry is simply 

_ 2_ 

'^=0 ^ ^ e e=0 

for points on the right side of the tunnel axis, or 


^ l£=o 


1 ^A ! 
p d 0 g=n 


for points on the left side of the tunnel exis. 


Inspection of equation (31) 
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shows Immediately that the 0-derivative of the second term in the "bracket 
is zero for either case and the contribution of the first term is 


l^A 

p 80 


1 

p 80 


Furthermore , all vertical velocities on the axis itself may be obtained 
by considering only m. = 1, because 


0=n 


- - \ ZZ ZZ wv\(i,p) 


m = 1 n = 0 


p-n 


o , z_ “V' 

r m - 1 TI — n 




1±m o P mn(^P) = 0 (m > 1) 

p — > 0 p 

The usual geometric symmetries also contribute toward simplifying the 
calculations. For example, if the horseshoe vortex lies in the hori- 
zontal plane of symmetry of the tunnel, 

= h^^) = 0 

If, in addition, the vertical plane of symmetry of the tunnel is also 
the vertical plane of symmetry of the horseshoe vortex, all even values 
of m are el imin ated; in the corresponding antisymmetrical case (as 
with aileron deflection) all odd values of m are eliminated. 


Cylindrically Symmetric Term 

For a normal velocity at the tunnel wall go(l) that is independent 
of 0 the potential function cannot be given exactly in the form of the 
preceding section since for m = 0 the integral with respect to q will, 
in general, not converge for q in the neighborhood of zero. It is 
necessary to add additional terms to the potential so as to insure the 
convergence of the integrals with respect to q. Moreover, these terms 
must be of such a nature that the potential function is still harmonic 
and gives the required normal velocity at the tunnel wall. 
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The singularity-free potential inside the tunnel then takes the form 


0 


i. 


0 


Jo(iPQ) /' 
iqJo'Uq) J_ ( 


go ( p) 008 q(P ■ 1) d P 


2rtk£ 

q 


2 _ 

q 2 





-00 


8q(P) d P 


dq 


where 


k = 


13m 


Sq(P) 


|Pl- 


value of the integral, 
to exist* 


and where 


8q(P) d P 


is the Cauchy principal 


Both the limit and the integral must he assumed 


The appearance of these additional terms is not wholly due to the 
presence of the open section in the tunnel. For a source in a completely 
closed tunnel the second term does not vanish and would have to he used 
in calculating the tunnel- induced perturbation velocity hy the method of 
reference 13* However, for a closed body or a vortex system plus its 
reflections in a completely closed tunnel, both of these additional 
terms vanish. 


It is easy to verify the fact that the additional terms do insure 
the convergence of the integral with respect to q. A straightforward 
differentiation then shows that $ is in fact harmonic and satisfies 


the boundary condition 



- 80 ( ^ ) • 


For the closed- open- closed tunnel, the boundary condition (see 
part I) that the velocities far upstream and downstream be equal is no 
longer automatically satisfied by putting the total tangential velocity 
on the jet surface equal to zero. The determining conditions for 0 A 

are now 


(1) 

(2) 

ao a 

8<d a 

dp 

II 

0 V 

11 a 

0 


(p < 1) 

( | < a and £ > b) 

(3) 

d$A 


d($Q + ®C) 

+ u 

P=1 

(a ^ £ < b) 

d£ 

11 

rH 

II 

a 

dF~ 
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00 

(5) 


lim 

+00 


w 


= lim 




= 0 


(I = a) 


Conditions (l) and (2) are satisfied by assuming go(0 = 0 for I 
and i > b . Thus 


0 


A 


1 

rt 



J 0 (ipq) 

iqJ 0 '(iq) 



g Q (p) cos q(p - 



Again it is desired to find go(t) for a < | < b so that 
satisfy conditions (3), (4), and (5 )* The representation of go(0 
the same form as before (equation 27) 


00 


goU) = h 00 sin | £-^-f + XZ kon sin m b^I 

n = 1 


automatically satisfies condition (5) • Substituting this series in 
equation (32) gives 


00 

^A = Z — h 0n P 0n^' p ^ 

n = 0 


where 


p 00 (tp) 


Jo L - 


JoUpq) 


W^M)J a sin 2b 


sin o f" : "! 008 q(P - 0 <ip 




dq 


< a 


(32) 


will 

in 


« 
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/ CO 

i 

*- 


iqJo ' (fq)/ a 3ln HrH 008 9(P - 5 > 46 


/ sin nj tP ~ — dp 

q 2 4 15 - a _ 


dq 


The resulting infinite series for the P()n(t,P) obtained by means 
of contour integrations are given in appendix A. 

Condition (3) then becomes 




dPon(^l) S (°o + <1> C) 

Tl = - hi 


+ u (a < ^ < b) 


p=l 


8P0n 8Po n 

But (see appendix A) lim — = - lim 

!-=► +u° - 00 

and so lim ^- 7 — = - lim t so that condition (4) becomes 

05 t-^_oo o& 


£ +°° 


d<J A 

lim h = 0. Thus it is necessary that 

|-^+oo 


h 00 , »0n & - (- 1 ) 11 ] 

rt / ' nit 

2(b - a) b _ a 


There exists a unique value of u for which the coefficients hQ n will 
satisfy this equation and it can be found as follows: Let 


where 


khi = ^On' + ^n' ' 


Z 


n = 0 


. , ^On(lA) _ ^0 + 0 C) 

n 0n ~tl “ 




p=l 
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and 


-i 

n = 0 a§ 


then if 


hpo 1 + uhoo* 1 + V (hQn* + uhQn") jl - (-l) n j 


= 0 


rt 

2(h - a) 


IUt 


n = 1 


b - a 


n = 1 


u = - 


00 V, I 1 


2 W + Z & - (-D n ] 

n = 1 


The coefficients h^ 1 and ho n * 1 are found by solution of sets 
of simultaneous linear equations, as described in the previous section. 

The function $ A ,! =ZI hQ n 1 *Pon(^p) is the perturbation potential 
which, when added to that of a uniform flow, gives the potential of the 
disturbance-free expanding tunnel described in part I and indicated in 

C$A 1 ! 

figure 2. The corresponding perturbation velocities — - — have equal 
and opposite values at <» and -oo. 


EXAMPLE 


As a somewhat simplified illustration, the problem of a semi- 
infinite unit doublet distribution (degenerate horseshoe vortex) along 
the tunnel axis was considered. The tunnel was assumed to have an open 
jet, 3 tunnel radii in length. The tunnel interference was calculated 
for four different positions of the upstream end of the doublet distri- 
bution, these positions being 0.1, 0.4, 0*7, and 1.0 radii downstream 
from the entrance . If the upstream end is taken as the origin of the 
coordinate system, then (see appendix C) 

r l <l)(S) * • (1 + 5 2)37 5 ' 1 1 0 iVUq)[?0 (g) + Kl(,3 3 oos 91 





(m ^ 1, J i 1) 


(I) = 0 


(33) 
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where Kq and Ky are the Bessel functions as defined in reference 15* 

The points JjjJ were taken as a, a + 0*3, a + 0.6, . . . a + 2.7, 
that is, a set of ten points, at 0.3 intervals, starting with the entrance 
lip of the tunnel. The coefficients h i n ^ were found ty the method of 


least squares for N = 0, 1, 


equations ri'-*-)(li) =2 h i n 


( 1 ) 


5 and 'also so as to satisfy the 

dPlnU,l) 


81 


at all ten points. Plots 


of the resulting functions g-^-^ ( |) for the different values of N 

indicated that convergence was essentially complete for N "between 3 
and 5 • This simplification results in appreciable saving in the amount 
of computation. Not only is it necessary to solve a smaller set of 
simultaneous equations, but also P ln and 8P i n /6^ need be found for 
fewer values of n. 


The computation was fairly straightforward. In the determination 
of r^U), Kq ^ K 1 were obtained from the tables of reference 15, 

and and J 1 ' from the tables and from the relations between the 

Bessel functions and their derivatives (references 15 and 16) . Weddle's 
formula (reference 17) is convenient for performing the integrations. 

In the case of P in , the values of ysl appearing in the formula 

for Qin^ s Hp) were found from the formula in appendix III of refer- 
ence 15, and J]_ and J-^ 1 as just noted. In the evaluation 

of lim — P 1ti it is noted that the value of lim — J-i (x) = 

_ O - L1J „ X M c. 

p — > 0 x — > 0 


The results of these computations, together with those for the 
completely open and completely closed tunnels and those given in refer- 
ence 3 are shown in figures 24, 25, and 26. In figure 24, the vertical 
tunnel- induced velocity along the axis for the four different positions 
of the lifting element together with the results for the open and closed 
tunnels are plotted against distance from the lifting element. The 
same results are plotted against the longitudinal distance from the 
entrance lip in figure 25 * Figure 26 shows the results of reference 3 
compared with the results of this paper for the same case - that of the 
lifting element 1 radius downstream from the entrance lip. 

The tunnel- induced velocity in the upstream regions and in the 
neighborhood of the lifting element, although only slightly less than 
that for an open tunnel for the lifting element 1 radius downstream, 
falls off more and more rapidly as the lifting element is moved towards 
the entrance lip. The maximum induced velocity is attained about 1 radius 
upstream from the exit, and is never more than j8 percent of the maximum 
value for a completely open tunnel. After the maximum the values fall 
rapidly and approach the values for a closed ‘tunnel in the downstream 
regions. The results of reference 3 (see fig. 26) are consistently below 
the present results especially in the region behind the lifting element. 
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An extrapolation from the present results indicates that the induced 
upflow at the lifting element, for the lifting element in the plane of 
the entrance lip, is approximately zero, or the average between the 
completely open and completely closed cases. The same result (that the 
effect in the plane of the entrance lip is the average of the effects 
for the completely open and the completely closed tunnels) was also 
obtained for the two-dimensional tunnel (fig. 18 ) . 


CONCLUSION 


For an open wind tunnel, the corrections corresponding to an 
infinitely long open jet will usually be adequately accurate if the 
region of interest (where the lift is located and where the boundary- 
induced flow is being considered) is at least half the jet height from 
the jet entrance and exit sections. As the distance of the lifting 
element from the entrance is decreased below this limit, the boundary- 
induced flow decreases rapidly and, when the lifting element is in the 
entrance plane, the induced angle at the lifting element is about the 
average of that for an open tunnel and that for a closed tunnel. 

In the theoretical studies of these flows, the usual boundary condi- 
tions of pressure uniformity on the free surface and of zero normal 
velocity on the closed surface must be supplemented with the conditions 
that the velocity be continuous at the entrance lip and that the velocities 
far downstream and far upstream in the closed sections be equal. For 
the two-dimensional open tunnel, a convenient general mathematical approach 
is to transform the infinite strip (representing the tunnel) to the upper 
half- plane by the logarithmic transformation and then to develop the 
desired complex velocity in this transformed plane. For the circular 
open tunnel the solution may be effected by expressing the potential by 
a finite series of Bessel functions, satisfying the boundary condition 
on the free surface at a finite number of points, and solving for the 
coefficients by simultaneous linear equations. 

For noncircular open wind tunnels, solutions in terms of available 
functions will be very inconvenient. For such cases, the trends indicated 
by the present results may suffice, when applied to the presumably known 
corrections for the infinitely long open and closed configurations, to 
provide adequate corrections. Solutions for the general three-dimensional 
configuration may also be possible by electrical- analogy methods, in 

which either the perturbation velocity potential or the acceleration 
potential is analogous to the electrical potential in an electrolyte 
solution. Such analogies may be characterized, however, by considerable 
technical difficulty. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., December 20, 19^3 
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APPENDIX A 

EVALUATION OF P^C^p) 

Evaluation for m ^ 0 

Contour integration .- For n ^ 0, 

n oo j (i DQ ) pi} 

= rt J i q j m i( lq ) d( 3 J 8in m |“E^ 008 q(P - 5) d P (Al) 

The inner integral may he found directly: 


sin nj cj* ~ - cos q(p - |) dP 
D ™ Q» 


1 

2 


p 

/■* 

p- _ 


r o ?t 

1 

sin 

a + q(P - 0 

+ sin 

“fh - i(p - 8) 

J a 

_ 



J J 


dp 


rut 


h - a 


f -(f^) L 


(-l) n cos q(b - - cos q(£ - a) 


(A2) 


The problem of evaluating P mn ( |,p) thus reduces to that of evaluating 
integrals of the form 

r J m (lpq) COS kq dq 
^0 iqJm'( i< 3)(^ 2 " h2 ) 

Consider the integral in the complex z-plane 

J m (ipz)e ikz 




izJ m '(iz)(z2 - h2) 


dz 
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* 


around the contour indicated in figure 27 • Its value is the sum of the 
residues of the integrand at its poles inside the contour. These poles 
are the values of z for which J m 1 ( iz) =0. The zeros of J m ' will 
"be designated y em j they are real and may he obtained from the formula 

in appendix III of reference 15 • The poles of the integrand then occur 

at iz = y sm , that is, at z = -iy sm » Since only the poles within the 
• . 
contour are desired, only the negative zeros of J m are considered. 

The residue of j ' |i ~ z ) z = ~ iy sm is > 


lim 

z > " iy sm 


z + ^sm 

V(iz) 


lim 

z — > -Iv™ 
J sm 


__ z + ^sm 

Jm'Uz) - VCysJ 


since Jm'Cygm) = 0$ hy the definition of the derivative this expression 
reduces to 


- i 

* * (y sm) 

The residue of the integrand at z = -iy S m is thus 


Jm(pysm) e^ 3111 
ysmXysm^ + ^)Jm' '(ysm) 


But the Bessel functions satisfy the relation 


V ♦ iV * (l - |)jm - 0 


so that at x = yam> where J m ' = 0, 


J m"(ysm) = ^^J m (ysm) 

ysm 


whence, finally. 


1 

2iti 


JmCipzJe^ 2 



dz . 



JmCpysmjs^^ysm 


+ h 


- yarn ) J m(ysm) 
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where the number of terms in the summation depends on the radius of the 
outer semicircle. 

For k ^ 0, if the radius of the outer semicircle is allowed to 
approach infinity in discrete steps so as to avoid the poles of the 
integrand, the integral over the semicircle approaches zero. The 
limiting values of the integrals along the two inner semicircles, as 
their radii approach zero, are readily determined ty the ‘usual process as 


i J m (-iph)e-^ 

hh 2 J m '(-ih) 


and 


i J m (iph)e lkl1 

4h2j m '(ih) 

These two terms may he combined, after noting that reversing the sign 
of the argument in J m and J m ’ merely reverses the sign of their ratio, 
to 


J m (lph) 

2h 2 J m '(ih) 


sin kh 


Equating the total integral along the infinite contour to the sum of the 
residues thus gives 


1 , 0 0 J m (ipq)e lkq dq 

iqJ m ' (iq) ^q 2 - h 2 ) 


J m (iph) 

2h 2 J m '(ih) 


sin kh 



Jm(pysm) e’ty^ysm 
+ h2 )( m2 " ^sm^)^m(ysm) 


where the y sm terms are now defined as the positive zeros of Jm' 

instead of the negative zeros (if J m '(x) = 0, so also does J m '(-x)). 
By equating the imaginary part of the left-hand term to the right-hand 
side, which is a pure imaginary, there results, finally 
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J m (ipq) cos kg dq _ iti J m ( iph) 

0 iqJm , (i<3)(q 2 " h2) ' 2h2j m '(ih) 


sin kh 


+ n 


Jm(pysm)e _k;7sm y 


sm 


S = (A 


2 , 

ysm + 


* 2 X : 


111 “ ygm ) J m(ysm) 


Expressions for , n ^ 0, m / 0 .- In the preceding development 

it was assumed that k 0 , which was acceptable with regard to equa- 
tion (A2) in view of the fact that the cosine is an even function of the 
variable. This essentially nonnegative value of k must be retained, 
however, in the final expressions for Pmn (i »p ) t 


mn 


1 

2 






a v ®U 




rut 


(-l) n sin l b - l l™ 
b - a 


- sin 


t - 


a | rut 


b - a 


nit 


b - a 


Jm(pysm)ysm 


s = 0 

. " l^ _a |y sm 


H 

OJ_ 

\ 

1 — 

( w\S] 

7 

ro 

T L 

y smy 1 

U - a )_ 

i 

i 

i 

!_» 

J m^ 


l"b"^ Asm 


(A3) 


Now for a ^ I ^b 


(-l) n sin I V * 1 ™ 

b - a 


= (-l) n sin 


b - a -(I - a) 


b - a 


mtj 


i - a 
-sin r u t ^ - _ a 


= -sin njt- 


U - a| 
b - a 
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For I < a. 


(-l) n 3to l- V t i m - -(-l) n Bln 

l) - a 


b - a + (a - t) 


b - a 


■rut 


= sin rut: 


a - t 
b - a 


= sin rur 


b - a 


For t > b. 


t , \n . b - I rut 

(-D SU1 J__ 


(-l) n sln ^ - a - (U 


Id - a 


Z_a 1 


rut 


= sin rut: 


I - a 
b - a 


= sin rur 


j - a[ 
b - a 


The first term on the right-hand side of equation (A3) is thus equal to 
zero for I < a or I > b. The desired expressions for film ( £ t P ) > n ^ ® ) 
are therefore , for a ^ I <. b , 


= 


, S - a 

— : — — — r sin rutf - 

rut 1 /. njt \ b - a 




rut 


Id - a 


QQ — — 


s = 0 


and, for £ a or ?^b, 


P «,p) - - 

mn b - a 


s = 0 


V (e) <e) 


|l-a|y sm _ (_-L) n e - |b-|| ysm 
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where 


Smn 


(s) 


^m(Pysm)ysm 


/ \ 2 


ysm + (^-r^ 

£ 

S 

m 

% 

^ — ' 


(n 4 0) 


f 

The calculations for n = 0, which follow similar lines, are not 
given here. The final formulas are, for a I < h. 


*m0 ^ *>p) 


Jm ( lp 2 ? rhn -) 

in T i/ in \ 

2(1 - a) J - ( 2 (b - a ) J 


sin 2 I ~ a 
2 h - a 



n: e "(^“ a )ysm 

2(1 - a) 


for I a. 


finf ) ( £ # P ) 


s = 0 


(p) 


ypim e 


-(t>-£)ysm 



e -(a-l)y S m 


and, for I > h, 


WM = 


s = 0 


<W) (s> (p) 


r sm 


5“ ( £“l)ysm 


n 

2(1 - 


a) 


5 " ( - _a )ysm 


where 


9mo^(p) = 


Jm(pysm)ysm 
)2“ 


jysm + |2(h - a) " ysm 2 )Jm(ysm) 


Evaluation for m = 0 

The evaluation of £,p) for m = 0 proceeds essentially as 

before with the difference that the contour of integration must avoid 
the origin. For n = 0, the final formulas are for a < i < h. 


6o 
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p oo(£>p) = " 


2(1 - a)l 

2 i.-a „ J 0 | 

r ip ^ 
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r 4 Jt 
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OJ 
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2 (1 - a) 
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(s) 


and, for | > 1, 
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For n ^ 0, the corresponding formulas are, for a <. £ <1, 
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for 5 < a. 


P On(^p) = - |7a - |) - (-l) n (b - |j] 

b - a 


~ %n (8) e- (a ' l)ys0 - (-D» e ‘0>-t)7BoJ 

and, for I > b, 

P On(^p) = ^ [(£ - a) - (-l) n (| - 13)1 

t> - a ~ 


irh V (S) R (5 ' a)7s0 - (-l) n e-( 5 -»)yaO 
8 = 0 L_ 
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APPENDIX B 


PROOF THAT 


S,l) 

IS BOUNDED AS n — a ^ I ^ b 


Differentiating the formula from appendix A for Pmn(i,p) and 
putting p = 1 gives 


dP (|,1) 

mn 


^m \ % 


. x 1 ( j nit V 
00 

+ ™ \~ 
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cos rut: 
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The second term of the right-hand member is of order — for large n so 
that 


n 


11m 
n — »a> 


1 , 1 ) 

~ n — 


J J i_SS N , 

nj m V h - a) | - a 

= lim r cos mtr - 


n — >‘<x> 




h - a 


The cosine factor of this expression merely oscillates between 1 and -1. 
For the remaining factor, it is noted from the asymptotic expressions 
for the Bessel functions (reference 15, PP* 59 - 6l) that J m (it) is 

essentially of the form as t — >00 from which it can be readily 

V2rtt 

iJ^ ' (it) a PP roac k es unity as t — »oo. 


shown that the fraction 


NACA TN No „ 1826 


63 


APPENDIX C 

DERIVATION OE EQUATION (32) 


Equation (32) ms derived for use in calculating some of the results 
given in reference 13, "but it was not explicitly stated and discussed 
in that paper. Because certain steps in its development are nob obvious, 
the present outline of its derivation is given. Familiarity with refer- 
ence 13 will be assumed. 


Certain difficulties arise in the treatment of the doublet line 
directly j so the result is found by considering a horseshoe vortex of 
finite span and letting the span approach zero. Equation (6) of refer- 


ence 13 gives the formula for 



(where $2 is defined in refer- 


ence 13) corresponding to a horseshoe vortex of strength T and 
span a having one trailing vortex along the tunnel axis and the other 
to the right of the axis. The procedure for the doublet consists of 
letting the yaw angle i be zero, expanding the radicals in ascending 
powers of a, and proceeding to the next step in the analysis, where a 
will eventually be made to approach zero. In the expansions, powers 
of o higher than the first may be neglected except where a occurs 
in the product la, since I takes on infinite values j furthermore, 
since for the doublet the field should be symmetrical about the vertical 


plane of symmetry (0 = 
immediately eliminated 


tj), unsymmetrical factors, 
. The formula for 

Ul 


as a cos 0, 
is thus 


may be 



-JE. lim 
a —- >0 


- sin 91 


jo 

1 + | 2 a 2 




i sin e 


a(| 2 + sin 2 0) \ \J 1 + |2g2^ 


I sin G 

W~ + 


sin 2 0 


cos 2 0 


x/TVT* CTT^jj 


(Cl) 


According to the procedure of reference 13, it is necessary to make a 
Fourier analysis of the three terms in the braces and then to insert the 
Fourier coefficients in equation (8) of reference 13 . 
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In the first term in the traces, the egression in the brackets is 
the first and only Fourier coefficient. Changing 5 to p and inserting 
the expression into the inner integral of equation (8) of reference 13 
gives 


n 00 

lim / j ■ -,- l , - - \cos q(p - I) dp 

°~ >0 J-oo \ \l 1 + p 2 a 2 \|i + py 


which, if integrated by parts, reduces to 


lim 
a — > 0 


p 00 

~ m 

1 

a 

L 

_( i + ? 2 ) 3/2 ( 

1 + ( S 2 . 2 ) 3 / 2 


sin 


q(p - *) dp 


The contribution of the first term in the brackets is 


U' y.fy ^ (q) sinqs 

(See reference 15, P* 52-) The fact that the contribution of the second 
term in the brackets is zero follows immediately, upon performing the 
change of variable p = Pa, from the Fiemann-Lebesgue lemma (reference 11, 
p. 172). 

The third term in the braces of equation (Cl) is converted as 
follows : 


I sin 0 | 1 

cos^0 

1 sin 0 

+ sin^0 \/l + I 2 

(1 + 1 2 ^/ 2 . 

(l 7 i 2 )3/ 2 


so that again the first and only Fourier coefficient is given directly. 
Inserting it into the inner integral of equation (8) of reference 13 
and integrating by parts gives 


/ , CO 

00 (1 + 7F 2 


cos q(p - I) dP = 


00 


q sin q(P - £) dp 
\|l + P 2 


= -2 qKo ( q ) sin q| 
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The second tern in the braces of equation (Cl) is not a one- term 
Fourier series. Its n^ 1 Fourier coefficient is given by a constant 
times 


I sin 0 sin n0 / 1 \ 

^ 0 cr(^ 2 + sin 2 ©) y \/ 1 + 


Inserting this expression in the inner integral of equation (8) of refer- 
ence 13, and reversing the order of integration gives 


'2jt 


sin 0 sin n0 



Vi + 


cos q(|3 - £) 


dp d0 


After substitution of p = 3 < 7 , the limit of the inner integral becomes 


lim 
a — >0 



P 

p^ + c^sin 2 © 




dp 


Integration by parts and elimination of terms in cr 2 reduces this 
expression to 


lim 

o — >0 


1 

q 



2p 2 + 1 

p 2 (l + p£) 3 /2 


sin 




which is zero, by the Riemann-Lebesgue lemma. 


Finally, then, for the unit doublet 



$ _ 

2 


2 sin 0 


m 0 1 °° J i(iqp) r~ / v . - / — 1 

* J iqji'(iq) N K 0 (q) K l (q J ql dq 




W p=1 = ~ J Q "iJi'(iq) &b«> + K 1^3 008 & d< 3 


J-L(iq) 
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The potential $q of a unit doublet line along the axis is 


- 




+ i 


n 2 + 



whence 


_ sin 9 

W P-1 (f + l) 3 ^ 


The flow of the visual reflection vortices for wings of finite span 
reduces, as the span becomes arbitrarily small, to a uniform upflow in 
the finite section of the tunnel and therefore contributes nothing to 
the longitudinal velocity. 


The coefficient of sin 9 in - — * i V ^g L 

be the expression given in equation (33) • 


p=l 


is thus seen to 
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(a) Infinitely long open tunnel. 



(b) Open Jet between closed entrance exit regions. 

Figure 1.— Illustrations for discussion of surface perturbation velocity 

in open wind tunnels. 
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(a) Contracting jet. The pressure on the free surface exceeds 
berth the upstream and downstream pressure^ hut is very 
close to the upstream pressure. 





(h) Expanding jet. The pressure on the free surface is less 

than either the upstream or downstream pressure, hut is 
very close to the upstream pressure. 

Figure 2.- Contracting and expanding jets (2 or 3 dimensions). 
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(a) Straight exit. 



(h) Bell-mouth exit. 



(c) Enclosed space beneath the lower free surface (two- 
dimensional tunnel) . 


Figure 3*— Spillage from the lower lip of the exit. 
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Figure 4.— Two-dimensional Jet with different pressures on the two 

free surfaces. 
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(a) Deformation due to image system; no downwash 
at the wing itself. 




(b) Undeformed upstream flow; downwash at wing 
is half of that at infinity. 

Figure 5 .- Two-dimensional open tunnel of infinite length. 
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(a) Two-dimensional vortex, showing current lines between the 

two plates. 



* 


(b) Two-dimensional vortex in a perturbation field having a horizontal 

velocity component. 



(c) Three-dimensional element of lift. 



(d) Horseshoe vortex of finite span. 

Figure 6.— Velocity— potential analogies for two— and three-dimensional 

lifting elements. 
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(a) Vortex on the center line. 




r©~i 
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(h) Unsymmetrically located vortex. 


Figure 7.— Velocity— potential analogies for the two-dimensional 

closed— open tunnel. 
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(a) Vortex on the center line. 



(b) Unsymmetrically located vortex (incomplete representation). 




(c) Expanding or contracting jet. 

Figure 8.— Velocity— potential analogies for the two-dimensional closed- 

open— closed tunnel. 
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(a) Lifting element on the center line. 


« 




(b) Unsymmetrical location of the lifting element. 

Figure 9*— Acceleration— potential analogies for the two-dimensional 

closed— open tunnel. 
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(a) Lifting element on the center line. 



(b) Tunnel arrangement that corresponds to omitting the additional 

short strips. 



(c) Expanding or contracting Jet. 


© 




(d) Curving Jet. 

Figure 10.— Acceleration-potential analogies for the two-dimensional 

closed— open— closed tunnel. 
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(c) Open section represented by only 
two longitudinal strips. 


Figure 11.— Acceleration— potential analogy for three-dimensional 
closed— open tunnel and two approximate acceleration-potential 
analogies for three-dimensional closed— open-closed tunnels. 

The closed-open analogy may also be considered as an approxi- 
mate analogy for the closed-open-closed tunnel. 
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(a) Representation of 

lifting element. 
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(b) Representation of* open boundary in 
closed-open-closed tunnel. 


Figure 12.- Acceleration-potential analogies as the 
difference between two velocity-potential analogies 
slightly shifted relative to each other. 
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Figure 13 •— Physical and transformed, spaces for two-dimensional 
closed— open tunnel of unit height. 
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Figure 14.- Physical and transformed spaces for two- 

dimensional tunnel of unit height with one exit boundary. 
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Figure 15»— Physical and transformed spaces for symmetrical two- 
dimensional closed— openr-closed tunnel of unit height. 
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Figure 16 . «« Maps of the functions w-^z) and *3(2) 
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Figure 17*- Tunnel -induced angle on axis of two-dimensional 
closed-open tunnel, with vortex at several locations alon« 
axis . ® 
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Figure 18.- Tunnel-induced angle on axis of symmetrical two- 
dimensional closed-open-closed tunnel, with vortex at several 
locations along axis. Length of open section is 1*5 times 
tunnel height. 
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Figure 19 .- Tunnel- induced angle on axis of two-dimensional 
closed-open-closed tunnel having one exit lip, with vortex 
at two locations along axis. Length of lower free surface 
is 1,5 times tunnel height. 


MCA TN No. 1826 


'8 



Figure 20.- Tunne 1 - induced horizontal velocity on axis of 
two-dimensional closed-open-closed tunnel having one exit 
lip, with vortex at two locations along axis. Length of 
lower free surface is 1.5 times tunnel height. Ordinate is 

(Induced horizontal velocitv\ hv 

V L ) jr> 


(Induced horizontal velocity^! 2h 
V V • 
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Figure 21*- Comparison of tunnel -induced angles on axis for three 
types of two-dimensional tunnels. Length of open sections for 
the closed-open-closed tunnels is 1*5 times tunnel height# 
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Distance from entrance, 5, units of tunnel height 

Figure 22.- Tunnel -induced vertical velocity v on axis of 
symmetrical two-dimensional closed-open-closed tunnel 
having additional velocities of -u^ and u^ on the 
upper and lower free boundaries, respectively. Length of 
open section is 1,5 times tunnel height. 


NACA TN No. 182 6 


93 


£ 



F igure 23r Closed - open - closed tunnel , sh owing 
coordinate systems . 



.8 1.2 1.6 2.0 2.4 

Distance from lifting element f t units of tunnel radius 
Figure 24- Tunnel- induced velocity para merer along tunnel axis for several 

positions of the lifting element in a dosed -open-dosed circular tunnel. 
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Figure 2 7r Path of complex confour integration. 


